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ABSTRACT 


One and two point statistical properties are measured 
in the flow fields of a coflowing turbulent jet. Two 
different confining surfaces (one flat, one with large 
curvature) are placed adjacent to the lip of the circular 
nozzle, and the resultant effects on the flow field are 
determined. The one point quantities measured include 
mean velocities, turbulent intensities, velocity and 
concentration autocorrelations and power spectral den- 
sities, and intermittencies. From the autocorrelation 
curves, the Taylor microscale and the integral length 
scale are calculated. Two point quantities measured in- 
clude velocity and concentration space-time- correlations 
and pressure velocity correlations. From the velocity 
space-time correlations, iso-correlation contours are 
constructed along with the lines of maximum maximorum. 

These lines allow a picture of the flow pattern to be 
determined. The pressures monitored in the pressure- 
velocity correlations are measured both in the flow field 
and at the surface of the confining wall(s). 

The velocity measurements are made with a laser Doppler 
velocimeter in conjunction with a phase locked-loop pro- 
cessor. The intermittency profiles are determined using 
a light scattering technique developed by Shaughnessy and 
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Morton. A 1/8 inch condensor microphone is used as 
a pressure transducer. 

The Reynolds number of the flow at the exit plane of 
the nozzle is 22,600. Two ratios of the jet velocity at 
the exit plane to the speed of the outer tunnel are used 
in the investigation and they are = 5.1 and = 10.88. 



Morton . A 1/8 inch condensor microphone is used as 
a pressure transducer. 

The Reynolds number of the flow at the exit plane of 
the nozzle is 22,600. Two ratios of the jet velocity at 
the exit plane to the speed of the outer tunnel are used 
in the investigation and they are- = 5.1 and X. = 10.88. 
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Chapter 1 


INTRODUCTION TO TURBULENCE - CHAOS OR STRUCTURE- 

Section 1.0 Introduction 

Whether one observes the flows of water in rivers, 
lakes, and seas, examines the motions of gases in inter- 
stellar nebulae, or experiments with laboratory flows 
through pipes or nozzles, the conclusion is readily 
reached that turbulence is the most frequent state of 
fluidic behavior. Turbulent flows are described as flows 
in which velocity, pressure, and temperature are random 
functions of -both time and space. ^ Generally, the image 
that is brought to focus in one's mind is of complete 
chaos and disorder. Special cases of organized motion' 
existing along with the randomness and chaos of the more 

classically defined turbulence were known to exist such 

% 

a.s, for example, vortex shedding from a circular cylinder 

7 

up to a Reynolds number of 10 . But in these instances, 

the structures were regarded as being closely tied to 

their geometric origins and not characteristic of "well 

(2) 

developed" turbulence. 

It has really only been in the past few years with 
the advent of improved statistical techniques and more 
careful observations that researchers have tried to search 
for understanding of turbulent flows that are not quite so 
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chaotic as previously assumed. Conceptions in the past 
of the nature of turbulence were closely linked with the 
available experimental tools. In the forties, the hot wire 
technique for velocity measurements was sufficiently de- 
veloped to test the assumptions of various 'phenomenological 
theories more directly. A complete rejection of these 
theories was the result. Advances in statistical tech- 
niques in the sixties resulted in some support being given 
to the turbulent viscosity concept as an engineering tool. 
With the continuing development in hot wire anemometry, 
experiments concentrated on studying the spectral distri- 
bution of the turbulent kinetic energy with a particular 
emphasis on the fine scale of turbulence. At the same 

time, various rate terms in the turbulent kinetic energy 

(3) 

equation were investigated. 

Lately, many researches have focussed more attention 
on shear flow turbulence than the "simpler” but more 
academic problem of homogeneous turbulence. It was known 
that large structures are important in the development of 
turbulent shear flows. The earliest and most decisive 
attempts to define the form of large eddies were made by 
Townsend and his students. ^ ^ 

Kline and Runstaler ^ and Brown and Roshko^ re- 
ported observations that have had a significant effect on 
the direction of subsequent turbulent research. They used 
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very basic optical methods and described that in turbulent 
flows of simple geometry , there is some order in the mo-r 
tion with an observable chain of events reoccurring random- 
ly with a statistically definable mean. 

As more and more evidence is obtained, the concept of 
spatially coherent structures has become increasingly 
reinforced. Also, the evidence has implied that the 
temporal behavior of these structures has a strong bearing 
on the development of the whole turbulent flow. 

Section 1.1 The Mixing- Layer 

The simplest free shear flow that can be considered 
is the region between two parallel streams moving at 
different speeds. In the shear layer, the flow is usually 
in a turbulent state. The majority of the existing know- 
ledge concerning shear layers either is conceived with 
the growth of small disturbances near the origin of the- 

layer or with the fully turbulent section farther down- 

/ o \ (a) 

stream. Michalke and Stuart ' have reviewed the flow 
mechanisms at the start of the mixing layer. In the fully 
turbulent mixing layer, Wyganski and Fiedler divided 
their measurements into three categories: (1) conventional 

average measurements (2) time average measurements in the 
turbulent and non-turbulent zones and (3) ensemble average 
measurements conditioned to a specific location of the 
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interface. They divided the mixing region into two parts, 
one on the high velocity side and one on the low velocity 
side. The interfaces on each side of the mixing layer 
moved independently of each other. Finally, the tur- 
bulence, which is not subjected to a constant mean shear 
(i.e. the gradient oscillates depending on the location 
of the interface) is not homogeneous and in general non- 
isotropic. Roshko et al^ 1 "^ - ^ 12 ^ looked at the effect of 
varying parameters of the mixing layer such as the velocity 
ratio, the density ratio, and Mach numbers. Shadow pic- 
tures of the flow revealed the presence of well defined 
large structures which have the appearance of breaking 
waves-. It was also found that the mean properties of the 
flow were controlled by the large coherent structures and 
independent of the small scale turbulence. 

The significant gap between the results of the work 

pertaining to the linear instability of shear layers and 

the descriptions of turbulent mixing layers was the focus 

( 13 ) 

of Winant and Browand. Flow vxsualizatxon was used 

along with hot film measurements of conventional mean pro- 
perties and statistical measurements aimed at defining 
ensemble-averaged properties of the structures observed in 
the visual studies. Unstable waves were found to grow, 
and the' fluid was observed to roll up into discrete two 

4 

dimensional vortical structures. These turbulent vortices 
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interacted by rolling around each other, and a single 
vortex was formed. This pairing process was observed 
to’ occur repeatedly, and it controlled the growth of the 
mixing layer . 

A theoretical study of the rolling up process to 

(14) 

which Winant et al referred was presented by Moore. 

For the first time, the chaotic motion which was a feature 
of earlier studies was eliminated and the details of the 
outer portion of the spiral were determined. This was 
accomplished by using a point-vortex representation to 
numerically study the evaluation of an initially plane 
vortex sheet. The rate of rolling up was also calculated. 

Section 1.2 Structured Jet Turbulence 

As the initial development of two dimensional and 
circular jets consists of two annular mixing layers form- 
ing around the central potential core, one expects these 
layers to behave in a similar manner to that described in 
the previous section. 

The main physical features of a subsonic jet are- de- 
scribed by Davies. At the jet orifice there exists a 

thin, cylindrical vortex sheet which rapidly develops into 
a layer of highly sheared strong turbulence. This settles 
down to an annulus of strongly sheared, self-preserving, 
intense turbulence surrounding a core of potential flow 
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after approximately one nozzle diameter. The statistical and 
geometric characteristics of the -flow remain similar for 
another four or five diameters. A second adjustment region 
follows where the potential core has vanished and the mean 
shear decreases rapidly. 

Crow and Champagne, in their investigation of the 

noise producing region of a jet, studied the evolution of 
orderly flow with increasing Reynolds number. The in- 
stability of the jet changes from a sinusoidal surface 
ripple to a helix and then to a train of- axisymmetric 
waves as the Reynolds number is increased from 100 to 1000. 

At a Reynolds number equal to 10,000, two axisymmetric 
structures which can be detected are surface ripples on 
the jet column and a train of large scale vortex puffs. 

Crow and Champagne, found that when a periodic surging of a 
given amptitude and frequency was imposed on the flow at 
the jet exit which caused the formation of a fundamental 
wave, the preferred mode with a Strouhal number of .3 was the 
most dispersive wave on the jet column, and least capable 

of generating' a harmonic. 

(17) 

A. J. Yule described investigation of coherent 
motions in the first ten diameters of round jets. The 
turbulent f lpw was dominated by the interacting and co- 
alescing of ring vortices. The large organized structures 
in the turbulent region of the flowing fluid differ funda- 
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mentally from the laminar ring vortices at the exit plane. 
These structures possessed a strong and well ordered 
three dimensional nature. 

Experimental and analytical' work aimed at under- 
standing the three dimensional interaction of discrete 
vortex rings and their stability characteristics have been 
presented by Kambe and Takeo v ' and Oshima. v ' The 
coalescense process and the observed and predicted growth 
of vortex core waves had basic similarities with occur- 
rences in continuously running round jets. 

The structure of jets which are subjected to a period- 
ic forcing has also been investigated by Girard, et. al.^*^ 
and Binder. ; Measurements indicated appreciable effects 
for forcing Strouhal numbers greater than 0.04. The two 
most important effects were a higher initial spreading 
rate and an associated periodic velocity component which 
decayed downstream. 

For further description of the information gathered 

concerning the flow of turbulent jets, please refer to 

C22) 

Catalano, Morton and Humphris. 

Section 1.3 Wall Jets 

While more and more attention is being focussed on 
the search for structure in turbulent jets and_ mixing 
layers, the flow fields of three 'dimensional turbulent wall 
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jets are virtually being ignored. Before progressing 
further, a definition of the term wall jet should be dis- 
cussed. The presence of both a solid boundary and a free 

boundary of the same fluid as the jet leads to the term 

% 

wall- jet which seems to have first been described by 
Glauert. 

As in a free or coflowing jet problem, in a wall jet 
there exist two main zones of the flow: (1) An initial 

zone of flow established immediately beyond the jet exit 
and (2) a zone of established flow (Figure 1) . The 
initial zone can be further subdivided into the potential 
core where the turbulence level is still relatively low, 
the developing boundary layer on the surface of the wall 
and the mixing or transition region which has been 
theorized to have a similar flow mechanism to that of a 
free jet. 

The zone of established flow can be further sub- 
divided into a region of free mixing where the transport 
processes are again similar to a free jet and the region 
near the wall termed the inner layer. The inner layer 
can be considered a boundary layer growing under a tur- 
bulent, decelerating, super-stream. The upper edge of 
the boundary layer is the vertical location of maximum 
mean velocity at any downstream position. 
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Figure 1 Schematic 


of a Wall Jet 


The two dimensional wall jet was first investigated 
( 24 ) ( 23 ) 

by Forthman in 1936. Glauert, ; Schwarz and Cosart, 
Sigalla, ^ 2 ^ and Myer et. al.^ 27 ^ added theoretical 
background and analysis to the problem. It was found that 
the mean velocity profiles, excluding the viscous sub- 
layer near the wall, were similar. Mura Hari v ' pro- 
posed analytical methods for predicting the flow develop- 
ment region and the bed shear distribution in a wall jet. 

(29) 

Dvorak determined an analytical method for turbulent 
boundary layers and wall jets in which the effects of 
large longitudinal surface curvature and the associated 
normal pressure gradients were included for the first time. 
The Reynolds stress terms in the equations of motion were 
approximated using an eddy viscosity approach based on 
the concept of intermittency . Blackwelder and Kovasznay 
obtained space time correlations with large streamwise 
separation in a turbulent boundary layer with a zero- 
pressure gradient. The use of conditional averaging tech- 
niques led to an average picture of the velocities in the 
interfacial bulges. 

The three dimensional turbulent wall jet was in- 

(31) 

vesti gated by Rajaratnam and Rani. Experimental ob- 

servations on the velocity and length scales, and wall 
shear stress for wall jets issuing from circular, square, 
triangular, and elliptic nozzles were presented." It was 
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found that the length scale in the transverse direction 

grew about five times as fast as the length scale in the 

(32) 

vertical direction. Narain analytically investigated 
three dimensional flows and found that the maximum velocity 
shows three regions of decay. The potential core region 
is followed by a region where the velocity decay is de- 
pendent on the shape and aspect ratio of the orifice. Fol- 
lowing this region is. a radial wall jet type maximum 
velocity decay region. Hill^^) looked at the effect of 
confining walls on a coflowing turbulent jet. It was 
shown that before the jet spreads to the wall, both the 
mean velocity and turbulent shear stress are self- 
preserving. Afterwards, only the shear stress is subject 
to this condition. 

From the above discussion, it is clear that the re- 
cent changes in thinking concerning the structure of tur- 
bulence have not been applied to any significant degree in 
the study of wall jets. In fact, very little information 
has been gathered concerning the turbulence velocity field 
of a 3-dimensional wall jet. Also-, even in, the measurement 
of mean velocity components nature of the flow was often 
times inferred analytically rather than actually measured. 

In view of the new approach taken in the measurement 
of turbulent jets and mixing layers, it is the purpose of 
this investigation to determine the effects of wall surfaces 
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on th'e flow field of an axisymmetric jet in a coflowing air 
stream. Of primary interest in this investigation is the 
comparison of both one and two point statistical properties 
measured at the same locations .in the flow fields relative 
to the exit plane of the nozzle, with and without the wall 
structures present. 

In order to best discern the effects of the wall 

structures, the velocity and concentration statistical 

properties of the unconfined turbulent jet must be well 

documented. This documentation was the subject of the 

. (34) 

author's Master's thesis. The investigation was 

limited to one Reynolds number (Re = 22,600). The statistical 
properties mapped included mean 'velocities, turbulent 
intensities, intermittencies, and autocorrelations and 
power spectra for both the fluctuating velocity and con- 
centration fields. The velocity measurements were taken 
with a laser Doppler velocimeter and the results compared 
to hot wire anemometer and pressure probe data where available. 
Mean velocity profiles as well as - the decay of the 

centerline .mean velocity were compared to the results of 
(35) 

Lawrence. Though there was an order of magnitude 

difference in the Reynolds numbers between the two in- 
vestigations, the agreement was quite good. Axial and 
radial distributions of turbulent intensities were compared 
to those measured by Davies and Lawrence. Agreement 
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again was quite reasonable. The conditionally sampled 
nature of the measurements led to the observance of some- 
what higher turbulent intensities than was previously re- 
ported in the outer regions of the flow. The intermittency 
factor profiles were compared to data presented by Sarni, 
Carmody, and Rouse. v 1 The agreement was quite good. 

A rise in the intermittency in the region adjacent to the 
lip of the nozzle indicated the presence of a vortex ring, 
at the exit plane. The appearance of peaks in several of 
the velocity power spectra was in agreement with Ko and 
Davies v and Bradshaw. ' ; The spectrum shape was found 

to be completely different in the potential core region 
as compared to the mixing region. In the potential core, 

the "bump" dr peak in the power spectra were explained by 
( 15 ) ( 38 ) 

Davies v and Bradshaw et al. ' as being caused by the 
passing of vortices downstream. In the mixing region, the 
high background turbulence level masked out the relatively 
small velocity fluctuations due to the vortex street. 

Whenever possible, comparisons are made with data 
presented in the literature. However, as stated previously, 
very little information concerning three-dimensional wall 
jets was found to exist;. 

One application of the research would be in the area 

( 39 ) 

of short take-off and landing aircraft technology. 

Presently, one design under consideration is termed Upper 
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Surface Blowing. The basic tenet of this approach is that 
the exhaust gases of the jet engine are turned by the 
presence of the aircraft wing's upper surface and any 
trailing edge flap. It has been found that a significant 
increase in lifting capability is realized, but the 
buffeting that the structure must endure is quite severe. 
Hence, there exists a need for more information for a 
configuration that can be classified a three dimensional 
turbulent wall jet. 
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Chapter 2 


LASER DOPPLER -VELOCIMETRY 

Section 2 . 0 LDV Background 

Since the first published account of the use. of the 
Laser Doppler Velocimeter (LDV) appeared in 1964 , much 
effort has been devoted to the LDV's development. The 
fact that classical probing techniques such as pitot 
probes, temperature probes, or hot wires suffer from 
problems including frequency response, probe interference, 
high temperature environments, low gas densities and poor 
spatial resolution has made the LDV seem very promising 
in its use as a diagnostic tool. Because of the many de- 
tails involved in the LDV as a system, the reader is re- 
ferred to the proceedings of several LDV symposia (41) 

as a beginning step. This section of the paper will try 
to incorporate some of the more recent published develop- 
ments and include some of the progress of ongoing work. 

Much attention has been, given to the analysis of the 

(42) 

reliability of LDV measurements. Kried developed an 
approximate technique for an estimation of errors in 
velocity measurements due to flow variation within the 
finite dimensions of the scattering volume. The technique 
is applied to both frequency tracking and individual 
realization LDV applications. In addition, the analysis 
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allows estimation of errors incurred when velocity mea- 
surements are made at locations sufficiently near the wall 
that part of the scattering volume is truncated by the 

wall. . A similar type of investigation was completed by 

( 43 ) 

Karpuk and Tiederman . They found that velocity var- 
iance in the data originates from three sources: velocity 

fluctuations caused by turbulence at. the center of the 
probe volume; spatial variations in the mean velocity a- 
cross the probe volume; and spatial variations in velocity 
fluctuations across the probe volume. Often times in the 
literature, the term "ambiguity noise" is used to describe 
the errors resultant from a finite sized control volume. 

New applications for Laser Doppler Velocimetry are 

(44) 

constantly appearing. Baker, Hutchinsen, and Whitelaw * 

used an LDV system to measure the instantaneous velocity 

within flames . Continued development in this area might 

allow the LDV technique to be used in the analysis of flow 

in the cylinder of an internal combustion engine. Wisler 
( 45 ) 

and Mossey developed a system capable of investigating 
the flow field within the rotating blade row of a low 
speed axial compressor. The compactness and flexibility 

A 

of the LDV technique was demonstrated by Munoz, Mocker, 

/ A \ 

and Koehler. ' Munoz, et. al . developed CO^ laser 
Doppler detection system which can be used to measure re- 
motely the true airspeed of an aircraft. Hughes and 
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(AJ) 

Pike considered several types of Laser Doppler Veloci- 

metry which are suitable in the remote measurement of 

wind velocity. The most significant conclusion reached 

was that CO laser systems have considerably greater 
2 

sensitivity than visible laser systems at ranges greater 
than a few meters. Design considerations for a 3-D LDV 
system for studying gravity waves in shallow water were 
reported by Hallermeier. The laser velocimeter system 

used three frequency-modulated light beams and one de- 
tector for measurement of the instantaneous velocity vector 
in reversing flow. 

Measurements in free and confined coaxial jets with 

(AS ) 

recirculation were made by. Owen * . Measured turbulent 

kinetic energy .and shear stress distributions indicated 
that nonequilibrium effects are dominant in the initial mixing 
regions and that higher order velocity moment turbulence models 
will probably be required to compute these flows . Velocity 
measurements of a slightly different sort were proposed 
and the results reported by Wang, Clifford, and Ochs.^^ 

Wind and ref ractive-turbulence. were sensed using crossed 
laser beams . The technique interprets the naturally oc- 
curring atmospheric modulation of the beam intensity 
(scintillation) to extract the average transverse speed 
and direction of air motion between source and receiver. 
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Recently, the LDV has been applied to a variety of 
* 

flows in both the transonic and supersonic regimes. These 

flows include jets, boundary layers, flow field around 

bodies, wakes, and cascades. Pfeiffer, vom Steim, and 
(51) 

Koch have reported measurements of mean velocities 

of a free jet issuing from a shock tube where the run 

time was of the order of one millisecond. The results 

indicated that the LDV can be used to measure accurately 

f52 ) 

transient phenomena in high speed flows. Yanta and Lee • 

(53) 

and Rose and Johnson have made boundary layer mea- 

surements at a nominal Mach number of three with no 

pressure gradients. The final application to be dis- 

(54) 

cussed here was reported by Flack and Thompson. In 

the report, measurements of velocities in the free stream 
(inviscid) part of a transonic flow are made. The effect 
of parameters such as signal-to-noise ratios, particle 
response, seeding rates, etc. were then documented. 

Section 2.1 Laser Doppler Velocimetry - Theory 

The general theory underlying the LDV technique is 
that the flow velocity at a point is determined by mea- 
suring the Doppler shift of monochromatic laser light 
scattered from small particles in the flow. To increase 
the amount of light scattered, the flow is often seeded 
by injection of particles of suitable size and density 
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which are capable of following the fluctuations in the 
velocity field. 

For the general scattering geometry, the Doppler 
shift is given by f Q = (kg - k L ) • Q 

where kg and k^ are the wave vectors of the incident and 
scattered radiation and Q = pi +• q] + rm 

represents the velocity Q in terms of components along 

> 

and perpendicular to the direction .of the mean flow, i. 

For the special case when the vector (kg r k^) is 
parallel to the x- direction then f D = (2P/\o)sin (B/2) 
where P is the angle between kg and k L and Ao is the wave- 
length of the monochromatic laser beam. 

The question for each application of the LDV is to 
determine what parameters affect its reliability and to 
determine the result of the biasing of the flow field 
information when it is found to be important. The biasing 
which is the error in the measurements made by the LDV 
technique, can be related to the basic parameters of the 
LDV setup and the properties of the velocity field, and 
the actual dynamics of the particles used to seed the flow 

In Appendix A, the effects of the finite sized laser 
(55) 

control volume on the measurement of mean and fluct- 
uating velocities are addressed. It is found that the 
summation over a finite volume yields a good estimate of 
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the mean velocity at the center of the probe under certain 

assumptions . The fluctuating velocities are seen to be 

influenced by any mean velocity and mean square fluctau- 

tions across the probe volume. The characteristics of the 

signal processing unit are also discussed as well as the 

quality of the Doppler signals detected. 

In Appendix B, the question of how closely the dioctyl 

phthalate, DOP, particles follow the actual flow field is 
(56 ) 

discussed. v 1 The particle lag of the DOP particles in 
the nozzle is found to be negligible. Secondly, the 
particles act as a low pass filter with the 3dB point of 
the rms response occuring at approximately 4x10 5 Hz. 
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Chapter 3 


EQUIPMENT AND PROCEDURES 
Section 3 . 0 Flow System 

A jet whose compressed air is marked with dioctyl 
phthalate (DOP) is mounted inside the test section of a 
low turbulence ievel subsonic wind tunnel. The nozzle 
with a contraction ratio of 14 to 1 over a length of 15.9 
cm has a Reynolds number of 22,600 and an exit plane dia- 
meter of 2.14 cm. The entire jet assembly consists of a 
diffusion section, a flow straightening section which 
includes two fine mesh screens and plastic coated soda 
straws, and the nozzle itself. 

With the mounting of the assembly inside the wind 
tunnel, the capability exists for varying one of the more 
important parameters which characterize the flow field - 
the ratio of the exit plane velocity of the jet to the 
uniform, parallel, constant outer flow velocity, ' X^ . 

For the majority of the results reported in this in- 
vestigation, the secondary flow in the tunnel was kept 
at 3.20 m/sec. This results in a value of 5.1 for X . . To 
ascertain the effects of changing the speed of the. outer 
flow, a second larger value of X^ was chosen (i.e. Xj = 10.88). 
It should be pointed out that the speed of the tunnel flow 
was constantly monitored by use of a pressure probe and an 
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electronic manometer (Barocel Type 1012) . 

' The 1 meter long test section of the wind tunnel has 
the dimensions of 215 mm by 305 mm. The walls of the 
test section are plexiglass with optical quality flat 
glass windows installed for better LDV performance. The 
floor of the test section was machined at several lo- 
cations which enables the jet assembly to be positioned 
anywhere in the range from approximately zero to ten 
diameters upstream from the fixed x location of the 
vertical laser "probe”. 

As was mentioned previously, the compressed air is 
marked with dioctyl phthalate (DOP) . The DOP in con- 
junction with a Royco Model 258 Smoke Generator produces 
a smoke with a mean particle size of 0.3 microns. The 
output of the smoke generator is mixed with clean air 
from the laboratory's compressed air lines. 

The wind tunnel, electronic equipment, and optical 
components, with the exception of the laser, are all 
placed inside a Lindgren double walled copper screen room. 
The copper screening filters out radio frequencies which 
are in the range monitored by the LDV equipment. 

Since the experimental design upstream of the jet 
is an important parameter, several methods were used in 
determining the actual conditions. First, smoke (DOP) 
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was injected into the outer tunnel flow and the region 
near the contraction section of the jet assembly was 
closely monitored. No evidence of separation of the out- 
er flow from the jet was observed. Secondly/ pitot-static 
tube profiles were taken for several x/2r locations. 

Two such profiles are shown in Figure 2 for X^ = 5.1. 

At x/2r- = 0/ the mean velocity is zero at the lip of the 
o 

nozzle and then reaches the free stream value at approxi- 
mately 1 . 5 diameters radially out from the lip . At 

x/2r = I, the velocity defect has decreased so that the 
o - 

smallest fluid velocity is about one half the value of 
the free stream velocity. This trend continues until 
at x/2r Q = 6, virtually no defect is apparent. 

Finally, a pressure probe, connected to 1/8 inch 
condenser type microphone was traversed along the out- 
side of the contraction section of the jet assembly. By 
the sound produced/ it is easy to distinguish between 
laminar and turbulent flow. Once again no evidence of 
the separation of the tunnel flow from the jet assembly 
was found. 

Section 3.1 Wall Surfaces 

Two different types of confining surfaces are used 
in this investigation. The leading edge of either of the 
wall surfaces is placed adjacent to the lip of the nozzle. 
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Care is used in the positioning of the walls so that a 
smooth transition is secured between the lower lip and 
the upper edge of the confining surface. The positioning 
is indicated in Figure 1. The first type of surface 
is a flat aluminum plate (henceforth denoted "plate") 

30.48 cm long, 30.48 cm wide, and .635 cm thick. Figure 
3 gives a top and side view of the plate, and supports. 

The leading edge is contoured so that it disturbs the 
outer tunnel flow as little as possible. The entire 
assembly fits snugly into the wind tunnel test section. 

In order to use a forward scattering mode of LDV operation 
when the plate is placed in the test section, a portion 
of the aluminum has been milled out and. replaced by an 
optical quality glass flat. This allows the scattered 
light to be detected below the surface of the plate. 
Actually, two flat plates were constructed. Since sur- 
face pressure-velocity correlations were sought in the 
investigation, pressure ports had to be mounted in the 
confining walls. The second flat plate was necessarily 
fabricated entirely out of aluminum with the pressure 
ports mounted flush with the supper surface. Figure 3 
also presents a schematic view of the pressure ports 
themselves as well as the pattern of their location in 
the plate. 
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The second type of confining surface used is a scaled 
down version of the wing used at NASA-Langley in their 
"Beach” test. The wing or curved wall surface (henceforth 
denoted "flap") is 15.24 cm long,- 17.78 cm wide and sweeps 
out an arc of 70 degrees. The flap is constructed from 
white pine, with several coats of stain and wax applied to 
smooth out the surface and make it comparable to the finish 
of the plate. Figure 4 shows the flap and its mount, which 
again are placed in the test section of the wind tunnel. 

Also, the locations of the surface pressure ports are shown. 

A second point to consider is that the flap does not 
extend across the entire test section. To determine what 
effects if any were the result of end effects, several 
measurements for the flow over the flap were made with 
extended flanges on both sides of the curved wall structure. 
With these flanges, the flap extended across the entire 
test section. No noticeable difference in the measurements 
was determined . 

The inclusion of the plate or flap in the wind tunnel 
test section necessitates the monitoring of their effects 
on the outer tunnel flow. A pitot static probe in conjunction 
with an electronic manometer was used to make mean velocity 
profiles, traversing from near the top of the test section 
down to the surface of the flap, then from the underside of 
the flap to the lower wall of the wind tunnel. The profiles 
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were taken at both values of and for both the flow over 

the plate and the flow over the flap. Examples of the re- 
sultant plots are shown in Figure 5. The outer tunnel flow 
was found to behave as a uniform stream outside the region 
affected by the jet flow field. Profiles of the outer 
tunnel flow with the flap in place but the jet off were 
also determined. Again, the tunnel flow remained uniform 
outside the nozzle and flap wake region. 

The coordinate system used in this investigation is 
shown in Figure 6(a). In Figure 6 (b) , the jet/flap 
assembly is shown mounted in the test section of the wind 
tunnel. The optical components will be discussed at a 
later time. 

Section 3 . 2 Two Color Laser Doppler Velocimeter 

The light source for both the velocity and intermittency 
measurements is a Spectra Physics Model 165 Argon Ion Laser. 
The diameter at the 1/e 2 points of the laser beam at its 
head is 1.5 mm while in the test section the diameter is 
approximately 2 mm. The noise level can be reduced to about 

-i' 

0.5% by tuning. Eight frequencies are present in the all- 
lines-operating mode with the two strongest (488.0 nm rated 
at 700 mw and 514.5 nm rated at 800 mw) being selected for 
use. In order to split the laser light into its different 
wavelengths, an extra dense flint glass prism is used 
(n = 1.655) . 


29 



o o 
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Several views of the different components of the two 
point measuring assembly are shown in Figures 7 (a) -(d). 

In the pictures, the plate is present in the test section. 
Since the plate is constructed by use of an optical quality . 
glass flat, a forward scattering (reference beam) mode of 
light collection can be utilized for the vertical beams. 

The resultant higher light intensity levels can be detected 
by using a Bell & Howell photodiode with an active surface 
of 1.27 cm. Forward scattering can also be used when the 
freely- expanding coflowing jet is investigated. However, 
when the model flap is in place, the receiving optics de- 
tector assembly must be mounted ‘on the top side of the test 
section and, hence, a back scattering mode is used. 

Since the intensity of light collected in the back 
scatter mode is significantly less than in forward scattering, 
a photomultiplier tube (EMI-95025) is used due to its greater 
gain potential. The spatial filter for the photomultiplier 
tube consisted of a 10 cm lens mounted immediately above • 
the upper test section wall and a 400y (diameter) pinhole 
located at the focal plane of the image formed. The choice 
of the 400u pinhole gives the largest signal to noise ratio 
for the Doppler signal. 

For the horizontally directed beams, a forward scattering 
(dual beam) mode can be employed for all three flow con- 
figurations. Once again, a Bell & Howell diode with an 
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Figure 7(a) Two-color LDV Setup 



Figure 7(b) Jet and Plate Installed in Wind Tunnel 
Test Section (Note outline of jet on 
surface of the plate) 
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Figure 7 (c) Dual Beam Receiver/Detector (Diode) 



Figure 7 (d) Back Scattering Receiver/Detector 
(Photomultiplier tube) 
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active surface diameter of 1.27 cm is used. To increase 
the signal to noise ratio of the Doppler signal, a spatial 
filter consisting of a lens with a focal length of 10 cm 
and an aperture with a 50 0y diameter located at the focal 
plane of the image formed is utilized. 

Both optical units or "boxes” are identical in con- 
struction and components. The optical components consist 
of a rhomb iod shaped prism, a beam splitter, a mirror and 
a lens. The rhomboid prism allows the rotation of the en- 
tire unit about the optical center. With this capability, 
the plane of the two intersecting beams can be adjusted so 
as to give transverse as well as longitudinal velocity 
information. The. beamsplitters are rated at 70-30 (i.e. 

70% of laser beam passes through while 30% is reflected) . 
The focal length of the lens is 30 cm which produces a 

control volume size of approximately .05 mm (diameter) x 1 
(55) 

mm (length) . For the reference beam mode of operation, 

a #2 neutral density filter is used to reduce the intensity 
of the beam which hits the photodiode surface. 

The entire optical assembly is mounted on an aluminium 
slide which provided movement of both laser "probes" in a 
horizontal direction perpendicular to the jet and tunnel 
flow axes. Additionally, the horizontal beam unit can be 
moved independently to allow a displacement of the control 
volume in all three directions. 
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The two color LDV arrangement proved to be extremely 
effective in determining the one and two point statistical 
properties of the flow. For one point properties, flexibility 
is inherent in the system because when an opaque surface is 
encountered with one "probe", the use of the other laser 
measuring device allows the experimentalist to always use 
a forward scattering technique. The data obtained from the 
different "probes" at the same location in the flow field 
were remarkably similar except in the cases when sharp 
gradients in the flow field properties occurred (Appendix A) . 

Two point velocity measurements require a very precise 
determination of the separation between the sensitive control 
volumes . Critical to this determination is knowing when the 
"probes" are measuring at the same point in the flow field. 

To accomplish this alignment, use is made of the Federal 
Scientific Ubiquitous Correlator. 

The fluctuating voltage (velocity) signals from both 
of the "probes" are correlated. A relative downstream dis- 
placement between the control volumes would manifest itself 
as a shift in -the maximum value of the cross-correlation 
from a zero delay time. Any lateral displacement would 
cause a skewness in the cross-correlation curve and/or a 
decrease in magnitude of the value of the covariance at zero 
delay time. This trial and error procedure, though somewhat 
tedious, allows very precise positioning of the laser control 
volumes . 
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Once the control volumes are precisely located, and 
the scattered light is collected by the photodiodes and/or 
photomultiplier tube, the resultant information is then 
analyzed by the phase locked-loop processor. The PLL 

is discussed in detail in Appendix C. The two PLL's used 
in this investigation are constructed of identical components 
and their respective output signals are in pha«?A w-i th each 
other . 

A schematic of the complete electronic equipment set-up 
is shown in Figure 8 . 

With the PLL properly set so that it is centered at 
the mean frequency of the flow, the signal from the detector 
is disconnected. A sine wave at the same frequency from a 
Tektronix Type 190 B Signal Generator is then connected to 
the PLL, and the mean frequency is read using a Hewlett 
Packard Type 5216A Electronic Counter. This procedure was 
adopted since it allowed the mean frequency to be determined 
much more accurately than was possible by other methos. Once 
the detector signal is reconnected, the voltage error signal 
from the PLL is utilized. Root mean square voltages are 
measured with a Disa Type 55D35 Rms Voltmeter. The voltmeter 
has a frequency range of 1 Hz to lOOKHz, and signals with 
crest factors as high as 15 can be measured. In this in- 

I 

vestigation, the compressed air jet is marked with OOP while 
the secondary tunnel flow is kept clean. The phase locked 
loop is triggered only if marker particles are present in 
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Figure- 8 Schematic of Electronics Setup 

















the control volume. Thus, due to the high Schmidt number 
(4x10 4 ) of the flow, -the phase locked loop determines the 
velocity field of the jet -only.-- Measurements.- taken in this 
manner are considered "conditionally sampled". 

To obtain autocorrelation and cross correlation 
curves and coefficients, the voltage error signal (s) are 
analyzed by use of Federal Scientific Model UC-202B 
Ubiguitous Correlator. From the autocorrelation curves, 
Taylor microscales and integral scales can be determined. 
The microscales are evaluated by- numerically fitting a 
parabola at the origin of the curve, while the integral 
scales are determined by a numerical integration of the 
area under the autocorrelation plot. A computer program 
was developed which enables the power spectral density 
function to be obtained from the autocorrelation curves 
by using a Fourier Transform technique. 

The Biomation Model 8100 Transient Recorder and the 
Nelson Ross Electronics Spectrum Analyzer were quite 
useful in determining- the quality of the Doppler signal. 

The transient recorder, which stores the digital equi- 
valent of electrical waveforms 'in a memory, was used to 
make histograms, to take pictures and to make plots of 
the individual Doppler bursts. 

For the intermittency measurements, the cumulative 
probability density mode of the correlator is used. The 
cumulative probability density (CPD)- is defined as the 
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fraction of time the signal is less than or equal to a 
certain voltage. If the CPD of the concentration signal 
was measured at zero voltage, one minus this fraction 
would be the intermittency. Actually, the voltage for 
which the CPD is 99% is measured in. order that only the 
percentage of time the concentration signal exceeds the 
noise is determined. 

Section 3.3 Pressure Measuring System and pu Correlations 
To determine the static pressure and the wall suf- 

( 57 ) 

face pressure, the system developed by Schroeder and 

Herling VJo; was used. The essential items included a 

1/2 inch condenser type microphone and a Nagra IV-SJ tape 

recorder. When the cross correlations were made between 

the pressure and velocity fluctuating fields, both signals 

were filtered (10Hz - 1000 KHz) before being processed 

in order to achieve a good signal to noise ratio. The 

filters used were DISA 55D25 Auxilary Units. 

A computer program was used that enabled the spectra 

of both the pressure and velocity to be obtained as well 

as the coherence between the two signals. Coherence is 

essentially the value of the cross-correlation coefficient 

as a function of frequency. If G- and G denote the 

11 12 

Fourier Transform of the autocorrelation function of 
Input 1 and Input 2 (in this investigation, pressure and 
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velocity), then the coherence, S 2 , is defined as follows: 

12 

i G ’ i 2 

6 2 * -I — 1-2J 

12 G G - 
11 22 

where G represents the 'Fourier Transform of the cross- 
12 

correlation function of Input 1 and Input 2. 

Section 3 . 4 Intermittency Measuring System 

Intermittency is customarily defined as the percent 
time the flow is turbulent at a given location in the field. 
In this investigation, the compressed air jet flow was 
marked with dioctyl phthalate (DOP) while the tunnel flow 
was kept clean. ' Since the Schmidt number, which is defined 
as the ratio of the kinematic viscosity of the flow to the 
coefficient of the diffusivity of the DOP smoke particles, 
is very high (10 4 ) the smoke remained confined to the jet. 
The only way that smoke particles would be at a location 
in the flow is if the turbulent jet is also present. 

The only way clean air could be marked by the DOP was 
by being entrained into the jet flow. Thus, outside the 
potential core region, the presence of smoke indicated the 
presence of turbulence. 

To determine the intermittency profiles, the system 

(59 ) 

developed by Shaughnessy is utilized. The design 

includes a lens-mirror assembly which permits the selection 
of virtually any size control volume from .1' mm by .1 mm 
by .1 mm to 1.5 mm by 1.5 mm by 1.5 mm. The lens has a 
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focal length of 20 cm. The detection system uses a PIN 
silicon photodiode with an active surface diameter of .5 
mm mounted behind a spatial filter. Thet spatial filter 
assembly is composed of a simple lens and an aperture 
located in the focal -plane of the image found. The lens 
of this spatial filter assembly is positioned for unit 
magnification . 

Section 3 . 5 Concentration Measuring System 

To make the concentration fluctuation measurements, 
a small modification is made to the two point LDV setup. 
The back scattering collector/receiver is kept the same 
but the current caused by the scattering of laser light 
by the particles is monitored. In the dual beam, for- 
ward scattering mode, the spatial filter is similar ex- 
cept for the addition of a 400u aperture. Also, the 
Bell and Howell photodiode is replaced by a PIN silicon 
photodiode with an active surface of .51 mm. It should 
be pointed out that -both receiver/detectors are positioned 
so that they only monitor the light being scattered by 
the DOP particles as the smoke passes through the control 
volume. Auto and cross correlation curves are 'obtained 
by use of the Federal Scientific Correlator. The power 
spectra are obtained by numerically Fourier Transforming 
the autocorrelations. 
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Chapter 4 


EXPERIMENTAL RESULTS AND DISCUSSION 

One and two point statistical properties are measured 
in the flow fields of a coflowing turbulent jet, a co- 
flowing turbulent jet and plate, and a coflowing turbulent 
jet and flap. The one point statistical properties in- 
clude mean velocities and turbulent intensities in all 
three directions, autocorrelations and power spectra of 
the fluctuating velocity field and the concentration 
field, and intermittency profiles of the passive ad- 
mixtive flow field. From the autocorrelation curves, the' 
Taylor microscale and the integral length scale are 
calculated using Taylor's hypotheses ^ . The mean 
velocity profiles for the flows over the flap and plate 
are compared to existing theories when possible. The 
effects of varying the ratio of the inner jet velocity to 
the outer tunnel speed are also noted. 

The' two point statistical properties measured in- 
clude space-time autocorrelations of the three different 
velocity fields, space time autocorrelations of the con- 
centration field and pressure-velocity correlations be- 
tween various ppints in the turbulent velocity flows and 
pressures measured in the flow fields at the upper surfaces 
of both the plate and the flap and at a sting mounted 
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static probe. The space-time autocorrelation coefficients 
are used to construct iso-correlation contours. These iso- 
correlation contours yield a considerable amount of in- 
formation concerning the effects of a confining surface 
on the structure in the near field of a turbulent jet. 

All velocity measurements are made .with a laser 
Doppler velocimeter, the intermittency data is -obtained 
by use of a laser light scattering technique, and the 
pressures are monitored with a 1/2 inch condensor type 
microphone . 

Section 4 . 0 Mean Velocities 

In Figure 9 , mean velocity profiles for the longi- 
tudinal component are presented for two downstream lo- 
cations (x/2r Q = 2,4), and at one vertical position 
(z/2r Q - 0.5, the coflowing jet centerline) for all three 
flow fields with A^ = 5.1. The ratio of the local excess 
velocity, U - U^g, to the excess core velocity, U Q - Upg/ 
is plotted versus lateral displacement (non-dimension- 
alized by twice the radius of the jet) from the center- 
line of the jet nozzle. 

Initially, consider the comparison of the width of 
the three flow configurations. The freely expanding jet 
and the flow over the plate are quite comparable in width 
both at two and four diameters downstream. However, the 
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flow over the flap has a significant increase in width 
compared to the previous configurations at both down- 
stream locations. Secondly, notice the effect of the 
confining walls on the decay of the core velocity. The 
presence of either the flap or plate serves to break up 
the mean velocity field in the potential core region more 
rapidly. The flap is, however, much more effective. Thus, 
as the turbulent jet flows over the flap, the velocity 
field is widened and decelerates much more rapidly. 

The observation concerning the effects of the two 

surfaces on the width of the jet is quantified in Figures 

10(a) and (b) . Here y^/ro is plotted versus x/2r Q for 

z/2r = 0.5. The quantity, y'., is defined as the lateral 
o 111 

distance from the centerline to the location where the 
mean velocity is the arithmetic average of its value on 
axis and in the secondary stream. Data for both ratios 
of inner jet velocity to tunnel speed (i.e. X^ = 5.1 and 
Xj = 10.88) are presented, and the results compared to a 
theoretical curve developed by Squire and Trouncer 
for the expansion of a coflowing jet (X^ = 5) . 

First, consider the findings for X^ = 5.1. While 
the agreement between the theoretical curve and experi- 
mental results for the' coflowing jet is quite good, one 
effect of the confining surfaces can be seen. For the 
flow over the flap, the value -y' increases much more 
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rapidly than is, the case for the other two configurations. 
With the plate in place , however, the value of ym is some- 
what smaller than that for the coflowing jet for the first 
several diameters downstream but eventually is quite 
comparable in magnitude. For = 10.88, the results are 

somewhat less clear. Though the value of y for the flow 

m 

over the flap still seems to increases the most rapidly, the 
rate of increase is considerably less than is the case for 
A j = 5.1. It is difficult to draw any stron conclusions, 
however, since the data is somewhat sketchy. Here, also 
the mixing widths of the unconfined jet and flow over the 
plate are quite comparable in magnitude. 

In Figures 11 (a) and (b) , the decay of the center- 
line mean velocity at z/2r Q = 0.5 for varying downstream 
locations is presented. The meaning of the data for the 
two jet/confining surface configurations should be 
discussed. As the flow exits the circular nozzle, the 
plate and the flap serve to transform the "top hat" 
mean velocity profile, characteristic of an ideal jet, 
into profiles more closely resembling wall jets. Thus, 
the velocity at the nozzle centerline would be expected 
to decrease at a different rate for the jet/plate and 
jet/flap cases as compared to the unconfined coflowing 
jet case. Secondly, considering the curvature of the 
flap, if this curvature were to cause a change in the 
direction of the flow, a decrease in the mean velocity 
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in the x direction would be the result. For the down- 
stream locations shown, however, the mean velocity has 
not changed its direction by any large amount from the 
x direction. For example, at x/ 2r^ = 4 with X^ = 5.1, 
total mean velocity vector acts at an angle of -7.1° with 
the longitudinal direction. 

For X j = 5.1, the effectiveness of the flap in de- 
celerating the flow in the x-direction is apparent. The 
more rapid decay of the flow over the plate as compared 
to the unconfined jet is also noted. When X^ = 10.88, 
the decay of the unconfined jet behaves quite similarly- 
to .the case when X ^ has the smaller value. The plate 
is slightly less effective in decelerating the flow while 
the effectiveness of the flap is reduced significantly. 

The two observations that have been made concerning 
the relative rates of decay of the centerline mean 
velocity and the increase in the mixing widths for the 
flows over* the flap with X^ = 5.1 and Xj - 10.88 are 
consistent and merit some discussion. The indication 
seems to be that, as far as the mean velocity field is 
concerned, as the value of X^ gets larger and approaches 
the value corresponding to a free jet (X ^ = °0 , the flap 
becomes less effective in both widening and decelerating 
the flow. 
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It should be mentioned that the location x/2r Q = 4, 
corresponds to the point where the flap starts to curve. 
Since the flow field, with A ^ = 5.1, increases in width 
much more rapidly, the energy required for the accelera- 
tion of the slower moving fluid particles would ne- 
cessarily cause a reduction in the momentum of the flow 
somewhere else. Hence, to some extent a widening of the 
flow field would be responsible for the more rapid de- 
celeration which has been observed for the case when 

A . = 5.1. 

3 

To explore the effects of the curvature of the flap 
for the flow with - 5.1, mean velocity components in 

. j 

all three directions are measured. These measurements 

are obtained for locations over both the flat and the 

curved segments of the flap. Mean velocity parallel to 

* 

the surface of the flap, U , versus perpendicular height 

above the surface of the flap, n, is presented in Figure 

12 . An interesting .observation that can be made con- 

★ 

cerning this series of profiles is that U /U has approxi 
ately the same maximum value at x/- 2r Q - = .4 as at x/2r Q = 2 
This seems to indicate a net transfer of momentum’ from 
the core of the jet down towards the surface of the flap. 
Also, for the first four diameters downstream, the flow 
seems to be tending toward the classical wall jet mean 
velocity profile. This will be discussed in more detail 
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Figure 12 Development of Mean Velocity Profile in 
Vertical Plane, y/2r Q = 0 

(Xj = 5.1) 
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in a later section of the report. However , after x/2r Q ~ 4, 
any similarity between the measured velocity profile and. 
the theoretically predicted self-preserving shape of the 
classical wall jet is no longer present. 

Mean- velocity profiles in the vertical/ W, and in the 
longitudinal U, directions for various downstream locations 
(x/2r Q = 5.5/ 6.0/ and 6.6) and at the vertical position/ 
z/2r Q = -0.57 are presented in Figure 13. The effective- 
ness of the flap in turning the flow is indicated by 
these profiles in both a qualitative and quantitative 
manner. At the centerline of the profiles/ for example, 
the ratio of U to W varies from 1.45 to 1.17 to 1.8 as 
the flow progresses downstream. These values of U/W 
would correspond to angles of 34.6°, 40.52° , and 29.0° 
with respect to the X direction. 

Figure 14 presents a vectorial diagram of the re- 
sultant mean velocities at the centerline (|y(/2r Q = 0) 
as the flow proceeds downstream. The angles noted are 
with reference to the x direction, and correspond to the 
direction of the mean velocity vector at each downstream 
location . 

Continuing with the discussion of the flow over the 
flap, a profile of the mean- velocity, W, in the vertical 
direction, z, for the downstream location x/2r. Q - 4 and 
z/2r Q = 1.85 is presented in Figure 15. Comparing this 
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Figure 13 Longitudinal and Vertical M* 
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Figure 14 Vectorial Diagram of Flow Over Flap 
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W profile with the one at z/2r Q = 0.5 yields the obser- 
vation that the jet seems to be rotating or rolling up 
as it spreads out over the surface of the flap. While 
at z/2r = 0.5, the maximum value of W/{U - is ob- 

O O i? o 

tained at the centerline of the profile. For z/2r Q - .185, 
this velocity ratio remains small and constant in the 
central region of the velocity field (i.e. |y|/2r o = 0) 
and reaches a maximum at approximately |y[/2r o = .90-. 

Hence, here an indication of the complexity of the three 
dimensional jet as it negotiates the curved wall surface 
clearly exists. 

Mean velocities in the lateral direction, V, are 
measured and are also presented in Figure 15 . Recall that 
the velocity measurements made in this report are con- 
ditionally sampled and thus the V velocities are somewhat 
larger than might be otherwise expected for unconditional 
measurements. The mean lateral velocities are determined 
for the unconfined coflowing jet, the flow over the plate 
and the flow over the flap for the same downstream lo- 
cation (x/2r Q = 4 and z/2r Q = 0.5) relative to the exit 
plane of the nozzle. While the maximum lateral velocities 
for both the unconfined flow and flow over the plate are 
comparable in magnitude, this is not the case for the 
flow over .the flap. In fact, the maximum value of the 
ratio V/(U Q ' - Up S ) has approximately doubled the maximum 
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value found for the other two flow configurations. 

• A comparison of the mean velocity profile for the 
flow over the plate and the flow over the flap to the 
theoretically predicted classical wall jet profile is 
presented in Figure 16. The velocity profiles in the 
center plane •( | y | /2r Q = 0) are tested for the fully de- 
veloped classical wall jet type similarity by plotting 
(U - U FS ) A u max -* U FS ) versus z/z%. Clearly for the 
data presented here, the "top hat" profile of the flow 
leaving the nozzle is still apparent. For any parti- 
cular profile, Umax and the length scale zh are obtained 
from the smooth curve drawn through the experimental 
points. The length scale, zh, is the vertical distance 
from the wall surface to the point where (U - U_ e ) = 0.5 

h 

(Umax - U ps ) and dU/dy is negative. 

The data presented in Figure 16 corresponds to three 
different downstream locations (x/2r o = 3.76, 4.0, and 
4.24) . Focussing initially on the mean velocity profile 
for the flow over the plate, it is clear that these plots 
are quite different in shape from the classical wall jet 
model. This is not surprising for it has been theorized 
that the wall jet actually is slower than the free jet 
in achieving similarity. Rajaratnam et. al. suggested 
that similarity will not be achieved until approximately 
ten to fifteen diameters downstream.. 
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Figure 17(b) Mean Velocity Profiles in Horizontal Planes, x/2r 
Camparison to Similarity Curve (X. = 5.1) 


Turning next to a discussion of the flow over the 
flap at x/2r Q = 3.76, the mean velocity profile is com- 
parable to the profile for the plate. However at x/2r Q = 
4.0, the mean velocities ploted are close to the theo- 
retical classical wall jet curve. After the flow has 
passed the actual curved portion of the flap (x/2r Q = 4.24), 
the closeness to theory is- no longer present. In con- 
trast to the flow over the plate, however, the mean velo- 
city profile does lose its "top hat" or flat characteristic 
much- more rapidly. Once again evidence is seen of the 
earlier breakup of the potential core region of the tur- 
bulent jet as it flows over the flap. 

To check similarity of the mean velocity profiles 
in the horizontal planes at any height, z/2r Q , from the 
wall, (D-Upgi/^Ui^^-Upgl ‘is plotted versus y/y in 
Figure 17. The quantity y ^y 2 is the length scale in 
the y or transverse direction and is equal to the lateral 
displacement from the centerline for which (U - U F g) - 


°- 5 < U m a x - D PS>- 

The mean velocity profiles in the horizontal planes 
for the unconfined coflowing jet, the flow over the flap 
and the flow over the plate at a given z location are 
compared to the curve predicted by Goertler v ’ for a free 
circular jet. The profiles for *all three configurations 
are closer to the theoretical curve in the horizontal 
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plane than is the case for classical wall jet comparison 
discussed previously. The mean velocity profile for the 
flap configuration is the closest to the similarity curve 
while the unconfined jet has the largest discrepancies 
for the vertical heights z/2r Q = 0.5 and z/2r Q = 0.375. 
However, closer to the surface it is the flow over the 
flap which deviates the most from Goertler*s curve. 

From this discussion, it is apparent that the flow 
over the plate tends toward the fully developed tur- 
bulent shape in horizontal planes faster than in the 
vertical plane. Also, the flap seems most effective in 
transforming the flat velocity profile characteristic of 
a circular nozzle. With respect to the unconfined co- 
flowing jet the mean velocity profiles are closer to the 
fully developed turbulent shape in the horizontal planes 

at z/2r = 0.125 and 0.25 than is the case at z/2r = 
o o 

0.375 and 0.50. 

Section 4 . 1 Turbulent Intensities 

Turbulent intensity is the ratio of the root mean 
square of the turbulent velocity fluctuations to the mean 
velocity at some reference location. In the investigation, 
the turbulence level is non-dimensionalized by the excess 
centerline velocity at. the exit plane of the nozzle (i.e. 

U q - — - Mathematically,' turbulent intensity in the x 

direction is defined as follows: 
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where 


urms = u 2 dt)^ 


It should also be noted that the turbulence inten- 
sities were corrected for ambiguity noise using the methods 
of Morton. The ambiguity noise was found to be 

approximately 1.8%. 

The distributions of turbulent intensities in the 
longitudinal direction and in the horizontal plane lo- 
cated at z/2r^ = 0.5 are plotted versus lateral dis- 
placement (i.e. |y|-/2r } from the centerline in Figure 

o 

18. This data is for - 5.1. For the case of the 
unconfined, coflowing jet, the • turbulent intensity in- 
creases from about 2.5% to 5% from x/2r ~ 2 to x/2r - 4. 

o o 

In contrast to this development is the more rapid in- 
crease in turbulent intensity when either the plate or 
flap are in place. At x/2r Q = 2, the turbulent intensities 
are all close in magnitude near the centerline. The 
largest value in fact, is maintained by the flow over the 
plate. At X/2r Q . = 4, however, the highest turbulent 
intensity is maintained by-the flow over the flap. In 
fact urms/(U^ - U FS ) is approximately five times the value 
it possessed at X /2r Q - 2. 
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Figure 18 Longitudinal Turbulent Intensities Profiles, z/2r Q = 0.5 
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Another important observation concerns the relative 
turbulent intensity value of the three flow configurations 
in the region 0.5 <_ jy[/2r Q <_ 1. Whereas the flap seems 

i 

to be amplifying the turbulence field in this region as 
well, the intensities for the plate flow are consistently 
lower than either of the other two cases. The presence 
of the confining plate seems to be diminishing the tur- 
bulent fluctuations occuring in the outer regions of the 
flow. 

Turbulent intensities in the lateral direction, 

(v - Ujg) , are shown, in Figures 19 and 20. In 

Figure 19, a comparison of the intensities for the differ- 
ent flow configurations at the same downstream location 
(x/2r Q = 4 and z/2r Q =0.5) is presented. Notice that 
the turbulent intensity in the lateral direction is much 
higher in the shearing layer of the unconfined jet 
than in either of the other two configurations. In 
Figure 20, (v~ 2 ) (U Q - U ps ) presented for the flow 
over the flap. It indicates the change in the turbulence 
level as the flow reaches the curved portion of the flap 
structure. 

The change of the ratios of the exit plane core 
velocity to the outer tunnel flow speed is seen to have a 
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significant effect on the decay of the centerline mean 
velocity. Next, the effect of varying A^ on the growth 
in the turbulent intensity at the centerline will be 
discussed. In Figure 21(a) and (b) , the turbulent in- 
tensity at y/2r Q * 0 and z/2r Q = 0.5 is plotted versus 

downstream location. x/2r . For both A . = 5.1 and 

o 3 

A„ = 10.88, the turbulence increases at about the same 
rate for both the unconfined coflowing jet and the flow 
over the plate. The magnitude is consistently higher - 
for the plate configuration. In fact, A ^ seems to have 
very little effect on the experimental data. 

For the flow over the flap, however, the value of 
the parameter, A ^ , is of considerable importance. In the 
case with A ^ = 10.88., the increase in the turbulence is 
at only a slighter faster rate than is the case for the 
other two configurations. However, with Aj = 5.1, the 
increase is much more rapid. This fact is quite con- 
sistent with the observations made concerning the decay 
of the mean velocity. At the location in the flow field 
where the jet is decelerated and widened at the most 
rapid rate,- the turbulence is also amplified greatly 
giving another strong indication that the flap serves to 
quickly break up the potential core flow. 
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Section 4.2 Longitudinal Autocorrelations 


If u = u(x,y,z,t) and u’ ■ u(x + Ax,yy z ,t) are the 
longitudinal fluctuating velocites at different downstream 
locations (x, x + Ax) but at the same lateral and vertical 
positions (y,z) in the flow, then, at a given time, the 
correlation coefficient, R uu ’ t is defined by: 


R 1 


uu 


-uu' 


(u 2 ) (u* 2 ) 


The autocorrelation coefficient, R, is given by: 


R = 


uu* 

h h 

(u 2 ) (u* 2 ) 


where u* = u(x,y,z,t +‘ At). Autocorrelations can be con- 
verted to longitudinal correlations by using Taylor ' s 
hypothesis/^ The basic tenet of this hypothesis is as 
follows: 

0 at- (B) 

where B is a statistical property. The assumption that 
is required is that the turbulent intensity not be too 
large, and that the flow be stationary. 

As was the case for the discussions of the mean velo 


71 




city and turbulent intensity, focus will be placed on 
both determining the effects of the two surfaces on the 
turbulent jet and also to .examine the importance of the 
value of the parameter, A ^ . 

Initially, recall the vortex street model for a tur- 
bulent jet. In this model, it is theorized that vortices 
are shed from the lip of the nozzle and move downstream. 

In the potential core region of the flow, the turbulence 
level is very low. Thus, it is possible to detect the 
relatively small effects of the fluctuating pressure 
field on the turbulent velocity. The resultant fluctuating 
velocities are due to the passing downstream of vortices 
being shed at the lip of the nozzle. These shedding 
vortices would give the autocorrelation curve a damped 
sinusoidal shape. The spacing of the peaks of the curve 
would correspond to the dominant "frequency in the flow 
field. 

Consider the autocorrelations of the axial velocity 
fluctuations taken at the locations |y|/2r Q — 0 and 
[y[/2r Q =.0.5 -and at the vertical height, z/2r o - 0.5. 

The autocorrelations are determined in all three flow 
field at both two and four diameters downstream of the 
exit plane of the nozzle. The lateral positions |y[./2r Q = 

0 and jy|/2r Q = 0.5 correspond to the centerline of the 
flow field and the lip of the nozzle. In Figure 22, 
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curves are presented for the unconfined, coflowing jet 
and the flow over the flap. These were taken at the 
centerline for both flow fields. Notice that for the un- 
confined jet, the autocorrelation coefficient has a 
damped sinusoidal behavior at both two and four dia- 
meters downstream. This is characteristic of the auto- 
correlation coefficient in the potential' core region of 
the jet. Also, there is a slight change in the delay 
time which corresponds to a change in the dominant 
frequency of the flow. 

With the flap placed in the flow, the sinusoidal 
nature of the autocorrelation curve is first reduced at' 
x/2r Q = 2 and then eliminated completely at x/2r Q = 4. 
This is consistent with the rapid increase which was 
noted previously in the turbulent intensities in what 
was the potential core region. As is the case in the 
shearing layer region of the flow, the increase in the 
turbulence level masks out the relatively smaller 
velocity fluctuations- which are caused by the passing of 
vortices . 

Still limiting the scope of the present discussion 
to the effects of the flap on the turbulent jet, auto- 
correlation coefficients are presented for the lateral 
location |y|/2r Q = 0.5 in Figure 23. Here, the measure- 
ments are made in the -mixing region ‘of both flow 
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Figure 23 Autocorrelations for Jet and Jet/Flap, |y|/2r 




conf igurations . Notice that there is little difference 
in the autocorrelations curves . ' The very high turbulence 
levels serve to completely mask out the velocity 
fluctuations due to the passing vortices and eliminate 
the damped sinusoidal nature of the autocorrelation 
curves measured in the potential core. 

Turning next -to a discussion of the effects of the 
plate on the turbulent jet, longitudinal autocorrelations 
are presented in Figure 24. At the centerline location, 
the autocorrelation looks quite similar to those deter- 
mined for the unconfined case. In fact, for the case of 
the flow over the plate, the peaks are slightly more 
pronounced though the dominant frequency is the same at 
two diameters downstream. At four diameter downstream, 
the increase in the turbulence level for the flow over 
the plate seems to be slightly more effective in masking 
out the pressure fluctuations though not nearly as 
effective as the flap. In the mixing region as has been 
discussed before, the autocorrelation curves are similar 
in shape for all three flow configurations. 

As is the case in other measurements reported, 
determining the effects of varying the outer tunnel flow 
speed on the various autocorrelation coefficient curves 
is also desired. To this end, the measurements of the 

s 

autocorrelations are presented for the case when Aj - 10.88. 
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Figure 24 Autocorrelations for Jet/Plate, z/2r 
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The format used for the following data is as follows: 
two sets of curves are shown for each of the three flow 
configurations, one set corresponding to the centerline 
location ( [y[./2r = 0), the second for the lateral lo- 

cation directly downstream from the lip of the nozzle 
(|y[/2r Q = 0.5). The data is presented in this somewhat 
different manner than was the case before in order that 
the change in the shape of the autocorrelation curves 
as the flow progresses downstream be more easily dis- 
cerned. 

In Figures 25 and, 26, autocorrelation coefficient 
curves for the unconfined, coflowing jet are presented. 

The characteristic damped sinusoidal nature is apparent 
at two and four diameters downstream at the centerline, 
though there is a slight shift in the dominant frequency 
and a small increase in the damping at x/2r Q = 4 . At 
eight diameters downstream, the autocorrelation curve has 
a typical shape for a fully developed turbulent flow. 

Out in the shearing region corresponding to the lip of 
the nozzle, the autocorrelation coefficients curves are 
devoid of any dominant frequency and remain essentially 
unchanged as the flow proceeds downstream. 

Turning next to the case of the turbulent jet flowing 
over the plate, the resultant autocorrelation coefficient 
curves are shown in Figures 27 and 28. Here, as is the 
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Figure 27 Autocorrelations for Jet/Plate, jy|/2r Q = 0, z/2r Q =0.5 


(Aj = 10.88) 



2 . 


Figure 28 


Autocorrelations for Jet/Plate, 1^1/2 


(X . = 10.88) 



case for the unconfined jet, the curves presented for 
the centerline location have the damped sinusoidal nature 
at both two and four diameters downstream of the nozzle. 
Also, there exists a slight shift in the dominant fre- 
quency and an increase in the. damping. ■ At eight dia- 
meters, the turbulent velocity field, again seems to have 
"lost” its structure. It is instructive to compare the 
autocorrelation curves measured for both values of X ^ . 

At the centerline, the dominant frequency in the auto- 
correlation coefficient curves is slightly greater when 

X. = 10.88, both at x/2r = 2 and x/2r = 4. Also, the 
j o o 

magnitude of the peaks in the curves are consistently 
larger at large delay times when X^ = 10.88. 

Next, consider -the comparison between the results 
for the flow over the plate to the unconfined coflowing 
jet when X^ = 10.88. The dominant frequencies are shifted 
to a slightly higher value and the magnitudes of the 
individual peaks are damped to a greater extent with 
the plate in the place. This is exactly the same ob- 
servations that are made when X ^ is equal to 5.1. In 
the shearing region downstream from the lip of nozzle, 
the autocorrelation coefficient curves are more typically 
turbulent _ in nature and are quite similar to the data ob- 
tained when Xj = 5.1. 
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Finally, autocorrelation coefficient curves are 
shown for the flow of the turbulent jet over the flap in 
Figure 29 and 30. Here, something quite interesting and 
maybe unexpected is occur ing. At two diameters down- 
stream from the exit plane, the curve again has a damped 
sinusoidal type behavior. However, x/2r Q = 4, the flow 
cannot be characterized by a dominant frequency and be- 
haves as if it is fully turbulent. This more rapid 
break-up of the potential core with the flap in place' 
was noted for the case when Aj = 5.1, but that similar- 
ity is somewhat surprising. Recall that the flow did 
not decelerate nor widen as quickly when Aj is set equal 
to 10.88. Also, the turbulent intensities at the center- 
line of the flow did not increase as- rapidly. Yet, the 
autocorrelation coefficient curves indicates that the 
turbulent structure of the jet has been similarly dif- 
fused in both cases. Thus, there exists a strong in- 
dication that the structure of the turbulent velocity 
field is influenced to a far greater extent than the 
change in the time averaged mean velocities would in- 
dicate. Downstream from the lip of the nozzle, the 
curves again assume their customary shape. Here, a 
similarity exists between the cases when A ^ equals 5.1 
and also 10.88. 
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Focussing on the comparison of the autocorrelation 
coefficient curves in the potential core for all three 
flow configurations when iU equals 5.1 and when A ^ equals 
10.88, though the shapes are quite similar, there is one 
difference which should be noted. At the larger delay 
times, the peaks of the curves are somewhat larger in 
magnitude . Thus , for the same flow configurations , the 
damping at large delay times seems to be a function of 
the inner jet to outer flow velocities. 

Section 4.3 Integral Scale s 

The longitudinal integral time scale in the x direction 
is defined in terms of the autocorrelation coefficient as 
follows : 

T — / fc *u(t) u(t') dt 1 
1 ' ° 

or 

T ~ f fc *R(t') dt’ 

1 o 

where t* is the time at which the coefficient first 
reaches the value of zero. To obtain the length scale, 
the integral time scale is multiplied by the local mean 
velocity, again using the "frozen" turbulence approxi- 
mation. ^ 

Thus L = U x T 

l l 
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Physically/ the integral length scale is approximate- 
ly the largest turbulent scale in the flow. 

The growth of the integral length scale for the un- 
confined jet for both values of A ^ is addressed in Figures 
31 (a) and 31 (b) . 

Figures 31(a) deals with the development of the 

integral scale at the lateral location corresponding to 

the lip of the nozzle while Figure 31(b) is concerned with 

the centerline growth. At the lip of the nozzle/ the 

integral length scale is found to grow linearly downstream 

. (3 3 ) 

as has been predicted by Ijaurence The scale for the 

unconfined jet with A^ = 5.1 grows somewhat faster than is 

the case when A. = 10.88. At the centerline of the flow 

3 

field, the integral scale grows similarly for both values 
of A j . This is not too surprising since the velocity of 
the outer tunnel flow would not be expected to have a 
significant effect on the turbulence in this region for 
the first several diameters downstream. 

For the flow over the flap, the development of the 
integral scale is quite different than is the case for 
the coflowing jet and is presented in Figure 32. It 
seems possible that the values measured at both the center- 
line of the flow field and out at the lateral location of 
the lip can be plotted on one curve. Also, the value of 
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X. does not seem to influence the results. The rapid in- 
3 

crease and then decrease of the values of the integral 
scales are notable in two respects. First, the large jump 
in the value of actually takes place before the flow 
has reached the curved • portion of the flap. This is also 
true with respect to the steep decline. An understanding 
of this behavior lies in an examination of the definition 
of the integral length scale. First, since the flow 
field over 'the flap grows at the most rapid pace, one 
would expect the integral length scale, which is an 
estimate of the largest scale in the flow to reflect this 
increase. However, as the flow approaches the actual 
curved portion of the flap, there occurs a turning of the 
flow direction and also a rapid transforming of the "top 
hat" mean velocity profile into a profile more typical 
of a wall jet. " Thus, the faster moving fluid particles 
are closer to the flap surface while out at z/2r Q = 0.5, 
the fluid is moving considerably slower. This large de- 
crease in the mean velocity would cause a decrease in the 
magnitude of the integral length scale if the integral 
time scale does not increase significantly since it is the 
product of integral time scale and the local mean velocity 
which are of importance. In fact, the integral time 
scale is nearly constant between x/2r Q = 2 and x/2r o = 4. 
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Secondly, as was noted previously, the integral scale 
is at best a very weak function of the value of A . 

» 

Though the presence of the flap amplifies the turbulence 
in the jet to a lesser extent when A_. is equal to 10.88, 
the product of the local mean velocity and the local 
integral time scale seems independent of the outer tunnel 
speed. This, however, is difficult to determine pre- 
cisely since the gradients of the integral scale are 
very large . 

The growth of the integral scales for the flow over 
the plate is presented in Figures 33(a) and 33(b). For 
both Xj = 5.1 and A ^ - 10.88, the integral scale at the 
centerline is consistently larger in magnitude than it 
is at the lip of the nozzle. After roughly four diameters 
downstream from the exit plane, the growth rate for all 
four curves shown is approximately the same. In, the 
data presented here, the value of A^ does effect the 
magnitude of the integral scale for the first few dia- 
meters downstream. However, after that, virtually no 
influence can be ascertained. 

Section 4.4 Taylor Microscales 

Let u 1 be the fluctuation of the longitudinal velocity 
in the ambient flow, then the Taylor microscale, A T , of 
turbulence is determined from the following condition. 
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Figure 33(b) Development of Integral Scale for Jet/ 
Plate, z/2r =0-5 (X. = 5.1) 



2 


( 



For isotropic turbulence/ the Taylor microscale can 
be rewritten: 

r 

1_ _ 1/5 J o k 2 S(k)dk 

* T 2 /” E(k)dk 


where k is a wavenumber, and E(k) is the three dimensional 
energy spectrum function- The integral in the numerator 
is an indicator of the overall scale of dissipation while 
the integral in the denominator is an indicator of the 
overall scale of energy. Thus X^ may be thought of as 
the ratio of the dissipation of energy to the total a- 
mount of energy in flow. 

The Taylor microscale can be determined experimental- 
ly by fitting a parabola to the autocorrelation coefficient 
curve near the origin. The intersection of the parabola 
with the delay time axis s along with the use of Taylor's 
hypothesis yield the' value of the Taylor microscale. 

Table 1 shows the value of the Taylor microscale for 
the three flow configurations with X ^ equal to 5.1 and 
also equal to 10.88. For each flow field, the Taylor 
microscale is determined both at the centerline and at 
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Comparison of Taylor Microscales 
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the lateral position corresponding to the lip of the 
nozzle. 

Consider initially the microscales determined for 

the unconfined turbulent jet. The values shown for. the 

centerline position are fairly constant though there does 

seem to be a slight tendency to increase. Downstream 

from the lip of the nozzle, that tendency is seemingly 

stronger. It is interesting to note that the values of 

A_ at the centerline when A. is equal to 5.1 are con- 
T 3 

sistently lower than correpsonding values when A ^ equals 
10.88. The opposite is the case, however, downstream 
from the lip of the nozzle. 

Next, the case of the flow of the jet .over the plate 
will be discussed. Here again, a slight increase in the 
value of the Taylor microscale is noted as the flow pro- 
gresses downstream. The value of A^ does seem to have a 
strong influence on the value of A^ downstream from the 
lip of the nozzle. With A^ equal to 10.88 the values of 
A t out in the shearing layer are almost an. order of 
magnitude different than for the corresponding position 
for the smaller value of A ^ . 

Finally, consider the flow of the turbulent jet over 
the flap. Here at the centerline and the shearing layer 
location, the Taylor microscale first increases and then 
actually decreases in value. This is the case for both 
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ratios of inner jet velocity to outer tunnel flow speed. 
Since the Taylor microscale may be thought of as the 
ratio of the dissipation of energy to the total amount in 
the flow, consider the implication of X T decreasing. A 
decrease in the magnitude of the Taylor microscale re- 
presents a decrease in relative magnitudes of the total 
energy in the flow -to the dissipation in the flow. Con- 
sidering the total energy aspect first, since the tur- 
bulence obtains its energy from the mean velocity, a de- 
crease in the mean velocity will eventually show up as 

a decrease in the total energy of the flow. With 
* 

respect to the relative increase in the dissipation of 
energy, recall that the dissipation is dominant at high 
frequencies. The flap which serves to amplify the tur- 
bulence in the flow field quite rapidly enhances the 
energy cascading process and, thus, could be . responsible 
for the relative increase in the magnitude of the 
dissipation. Evidence thus points to the fact that the 
flap has a very strong influcence on the turbulent 
velocity field. 

Section 4 . 5 Power Spectral Densities 

To determine the spectral density curves, a computer 
program is used to take the Fourier Transform of the 
autocorrelation functions. The power spectra are nor- 
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malized so that the area underneath each curve is 1. 

Initially, data is presented for the three different 
flow configurations with A set equal to 5.1. In Figure 
34, the spectral density curves both at the centerline 
( |y|/2r Q = 0) and corresponding to the. lateral position 
of the lip of nozzle (jy|/2r o = 0.5) .are shown for a down- 
stream displacement of two diameters and with z/2r Q .= 0.5. 
The power spectrum in the potential core of a turbulent 
jet is characterized by a large "bump" or peak. 0 &) 

At the centerline the peak of the unconfined jet is clear- 
ly evident. In fact, as would be expected from observing 
the autocorrelation curve for this location the power 
spectra of both the flows over the plate and flap are also 
characterized by a large peak. For the flow over the 
plate, the peak has approximately the same magnitude as 
the unconfined flow but with a slight shift in frequency. 
With respect to the flow over the flap, the peak has been 
markedly diminished and, again, shifted to an even higher 
frequency. In the shearing layer, the distinctions among 
the three configurations are no longer evident. In -this 
region, the peak is virtually entirely masked out by the 
relatively high turbulence . levels . 

Power spectra for the downstream location of x/2r o = 

4 and z/2r Q = 0.5 are reported for the three flow con- 
figurations at the centerline and nozzle lip lateral 
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positions in Figure 35. Fo.r the unconfined jet, the 
potential core region of the jet still exists and is 
characterized by a large peak, slightly attenuated and 
shifted to a higher frequency than was the case at x/2r Q = 
2. For the case of the flow over the plate, the peak is 
also present though with a smaller magnitude than for 
the unconfined jet and at a slightly higher -frequency. 
Turning to the flow over the flap, the strong effect 
the curved wall surface has on the turbulent velocity 
field is clearly demonstrated. Notice that the char- 
acteristic peak is totally absent. 

Finally just as was seen two diameters downstream 
of the lip, the power spectra in the shearing region 
are quite similar for all three flow field configuration. 

With the purposes of viewing the change in the £ower 
spectra for each configuration from a slightly different 
perspective and also examining the effects .of -varying the 
ratio of the inner jet velocity to the outer tunnel speed, 
power spectral densities for the three flow field are 
presented with Aj set equal to 10.88. The format chosen 
for display of the data is identical to that used for 
the autocorrelation coefficient curves with A^ equal to 
10.88. For each flow configuration, one set of power 
spectra corresponding to the centerline location (|y|./2r Q 
0) is shown as well as a second set for the lateral 


100 






location directly downstream from the nozzle lip ([y|/2r Q = 
0.5). 

Focussing initially on the development of the un- 
confined coflowing jet , the characteristic large peak in 
the spectrum is clearly evident at x/2r Q •= 2 (Figure 36) . 
This peak is slightly attenuated and broadened four dia- 
meters downstream of the exit plane. At x/2r Q = 8, the 
power spectrum has lost its peak and seems to be approach- 
ing a fully developed turbulent nature. Turning back to 
the data presented for the flow with A^ equal to 5.1, 
the peak is at the same frequency and the relative 
decay of the "bump'' in the spectra behaves similarly. 

In Figure 37, the change in the power spectra 
in the mixing region . is shown. There seems to 
be a shift of energy from the high frequency components 
to the lower frequency components of the spectrum. Here 
again this is the same behavior found when A.. = 5.1. 

The development of- the power spectra for the flow 
over the plate is shown in Figures 38 and 39. At the 
centerline location, the power spectrum had a peak similar 
to that found for the unconfined jet. However, as was 
noted for the * case when A^ = 5.1, the peak is both 
attenuated and shifted to a higher frequencies. Notice 
that again, the peak in the power spectrum broadens as 
the flow progresses from two to four diameters down- 
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Figure 39 Downstream Development of Power Spectral Density 
for Jet/Plate, (y|/2r Q = 0.5, z/2r Q = 0.5 

(X. * 10.88) 

D 
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stream from the exit plane, and at x/2r Q = 8, the "bump" 
is absent. The dominant frequencies in the spectrum when 
= 5.1 and when A = 10.88 are the same. Out in the 
mixing region (Figure 39) there also is a shift of power 
from higher frequency components to the lower frequency 
components of the spectrum. This behavior is quite 
similar to that exhibited by the unconfined jet. 

Finally, power spectra for the turbulent jet flowing 
over the flap are presented in Figure 40 and 41. As 
might be expected, the large peak in the power spectrum 
is absent at x/2r Q = 4 at the centerline of the flow 
field. However, there is one difference for the case 
where A_. is. equal to 10.88 and this occurs at x/2r Q = 2. 
With Aj equal to 5.1, the peak in the spectrum is slight- 
ly attenuated and the frequency is shifted to higher 
frequencies as compared to the case of the unconfined 
jet. For Aj equal to 10.88, neither of these changes 
occur. In the mixing region of the flow field, the 
shift of energy .from high to low frequencies is apparent. 

Section 4 , 6 Intermittency 

The percent time of the occurrence of turbulence at 
any point within the flow field is called the inter- 
mittency factor (I) . It has been demonstrated by Van 
Atta/ 63 ^ Shaughnessy ^ and Catalano, et. al. that 
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the intermittency of a given flow field can be determined 
by monitoring the diffusion of any seeding particles as 
long as those marking particles are passive and provided 
the Schmidt number is relatively high. The Schmidt 
number, which is approximately 10^ for this experiment, 
is defined as the ratio between the kinematic viscosity 
of the fluid and the diffusion coefficient of the 
seeding particles (in this case, DOP) . Thus the Schmidt 
number relates viscous diffusion to mass diffusion. 

Profiles for intermittency are presented for the 
region outside the lip of the nozzle. The radial and 
axial distributions of the intermittency factor at 
various vertical locations for the unconfined jet flow 
and the flows over the flap and plate are given in Figure 
42. All the data shown is for the case when A.^ is equal 
to 5*. 1. 

For the vertical location z/2r Q = 0.5, the flow 
field with the flap in place is significantly wider than 
either the unconfined jet or the plate case. It is 
interesting to note, however, that after the very large 
initial increase in the width, the flow field actually 
decreases in width as the turbulent jet progresses down- 
stream. Thus as the flow starts to be turned by the 
flap, the width of the flow field at z/2r Q = 0.5 decreases. 
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For the vertical position z/2r o =0.18, the turbulent jet 
flowing over the flap is greater in width at each down- 
stream location than it was at x/2r = 0.5 and it also 

. o 

increases in size as the jet progresses downstream. The 
flap seemingly acts to draw the jet down towards the 
surface and increase the flow field’s width. It also 
serves to transform what was an axisymmetric jet into a 
flow characterized by a rather elongated oval type 
cross section. 

In the discussion of the mixing widths for the flow 
fields, it was mentioned that they were about the same 
for both the unconfined jet and the flow over the plate. 
Here that discussion is reinforced by comparing the inter- 
mittency profiles at z/2r Q = 0.5. The profiles are quite 
similar in shape. 

Concentration Auto and Cross Correlations 

The concentration auto and cross correlations are 
defined in a similar manner as those determined for the 
velocity field. However, it is the fluctuations of the 
passive admixture (DOP) which is being monitored rather 
than the turbulent velocity fluctuations . The auto- 
correlation function of the concentration field is de- 
fined as follows: 

« _ ~cc*~ 

C (c 2 )' 2 i^ 2 ) 2 
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where c = c(x,y,z,t) and c* = c(x,y,z,t + At) 

Consider the cross-correlation function of the 
concentration field defined in the restricted manner 
such that only lateral or y displacements are allowed be- 
tween the two measuring points. Then the cross cor- 
relation takes the following form: 

c (x,y,z,t) c (x,y + Ay,z,t + t) 
l 2 

Rc c = t- r- 

1 2 (“ 2 ) % ( c” 2 ) 2 ■ 

1 2 

where c and c are the concentration fluctuations at two 
1 2 

different points in the flow field. 

Next, the further restriction is made that the 
concentration fluctuations will be monitored only at the 
lateral positions corresponding to the lip of nozzle 
(y/2r Q = ± 0.-5) at the vertical location z/2r Q = 0.5. 

Autocorrelations and cross correlations are measured 
in the three flow configurations at different downstream 
locations. Consider initially the case of the unconfined 
coflowing jet. Data is presented for x/2r Q - 4 and 
x/2r Q = 8 in Figures 43 and 44, respectively. At x/2r Q = 4, 
the autocorrelation of the concentration fluctuation is 
quite similar to the velocity spectia discussed earlier 
in this investigation. The curve, demonstrates a damped 
sinusoidal behavior. The cross correlation curve is 
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Cross-Correlation 


Rc c 
l 2 


-.05 








also of a similar shape. The maximum absolute value of 
the coefficient (Rc^c^ = .12) occurs at approximately 
zero delay time though there is evidence of a slight 
phase shift between the two concentration signals. At 
x/2r Q = 8 the autocorrelation curve is of a more typically 
fully turbulent nature. The maximum absolute value of 
the cross correlation coefficient has significantly de- 
creased (Rc c = .05) though there still is some relation 
1 2 ' 

between the fluctuating concentration field at y = ± 0.5. 

A large shift in the phase angle between the signals has 
also occured. 

Auto and cross correlations for the flow over the 
plate are presented in Figures 45 and .46 for the downstream 
positions, x/2r Q = 4 and 8. -The autocorrelation co- 
efficient curve at x/ 2r Q = 4 is quite similar to that 

found for the case of the unconfined jet at x/2r = 4. 

o 

The maximum absolute value of the cross correlation co- 

i 

efficient is similar to the previous case also (Rc c = .13) . 

l 2 

• The delay time for the maximum coefficient value has 
been shifted positively. The concentration signals are 
considerably more out of phase than is the case for the 
unconfined jet. At x/2r o = 8, the damped sinusoidal nature 
of the autocorrelation coefficient curve has disappeared. 

The maximum absolute value of the coefficient has de- 
creased (Rc c = .04). As in the case of the unconfined 
I 2 
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Figure 45 Auto and Cross Correlations for Jet/Plate Concentration Field 
x/ 2r o « 4, y/2r Q = ±0.5, z/2r Q « 0.5 ( A ^ = 5.1) 



jet, though the concentration autocorrelation curves 
appear to fully turbulent in nature, there still is 
evidence of a relationship {albeit a weak one) between 
the fluctuations at each lip of the nozzle. 

Finally, consider the flow of the jet over the flap. 
Auto and cross correlations at x/2r Q = 4 are presented 
in Figure 47. The autocorrelation coefficient is quite 
similar in shape to the previous two cases for this down- 
stream location. The maximum absolute value of the 
cross correlation coefficient is- also comparable 

(Rc c = .12). The cross correlation coefficient curve, 
l 2 

much like the case for the flow over the plate, demon- 
strates a large phase shift between the concentration 

fluctuations at y/2r = ± 0.5. 

o 

Concentration Field Spectra 

A spectrum of concentration or contaminant fluctuations 
is produced when smoke particles are mixed by the turbulent 
flow field of the jet. Many of the same concepts that 
are used in describing the kinetic energy spectrum can 
be applied to the spectra of the mixing smoke particles. 

The scales of the smoke particles fluctuations range from 
the scale of the energy-containing eddies to a smallest 
scale that is dependent on the ratio of the dif fusivities 
(Schmidt number) . The spectra are normalized for unit 
area. 
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Figure 47 Auto and Cross Correlations for Jet/Flap Concentration Field, 
x/2r Q a* 4, y/2r Q = ±0.5, z/2r Q =0.5 (Aj = 5.1) 


Concentration spectra measured for the three flow 

configurations are shown in Figures 48 and 49. In Figure 

48, the laser "probe" is positioned at x/2r Q =4, 

jy|/2r =0.5 and z/2r = 0.5. Consider first the case of 
o o 

the unconfined coflowing jet. Notice the relatively 
large peak in the spectra at approximately 300Hz . With 
the plate placed in the flow field, the peak is attenuated, 
and there also is an. increase in the energy at lower 
frequencies. By replacing the plate with the flap, this 
same trend is continued. The peak is diminished further 
and the energy at the lower frequencies' are slightly 
increased. 

In Figure 49, the concentration spectra for the jet 

and jet/plate flow are presented for the following 

downstream location: x/2r = 8, |y|/2r - 0.5 and z/2r = 

^ • o ' o 

0.5. At eight diameters downstream, the peak in the 
spectra has disappeared. Here again, the presence of the 
plate in the flow field has caused an increase in 
energy of the concentration spectra at low frequencies. 

Section 4.7 Space-Time Two Point Velocity Correlations 
The general space-time correlation coefficient for 
two points separated in a given flow field is as follows: 
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R , (x,y,z,t , % ,1 ,£ ,x) = 

11 1 ' 2 3 

u(x,y,z,t) u' (x+5,y+?,z + £?t +:t) 

1 2 3 

a* 

where 5 , t 5 '/and •£ represent the spatial displacements 
12 3 

of one point relative to the other and t represents the 
delay time. Here, the space- time correlations are limit 
ed in scope to fluctuating velocities in the x direction 
only. ' From these measurements, it is possible to cal- 
culate isocorrelation curves for optimum delay, 

Optimum delay time refers to the delay time at which 
for a given spatial separation the correlation is a 
maximum. 

Typical results obtained are shown in Figures 50 

thru 55. Figures 50 thru 52 deal with’ the case when 

5 is set equal to zero and hence are cross-sectional 
2 

views of the three flow fields in the x-z plane. Figure 
53 thru 55 correspond to the isocorrelations curves 
obtained with £ set equal to zero and thus represent 
cross-sectional views in the y-z plane. 

The reference point which remains fixed for all the 
correlation coefficients determined in the three flows 
has the following coordinates: 
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x/2r = 4 ! Y i /2r * 0 z/2r = 0.5 

o o 

This downstream location is chosen for several reasons. 
First, this is the position at which the flap first be- 
gins to curve in a concave manner away from the x-direction. 
Secondly, it has been observed in many of the one point 
measurements- presented earlier in this report that there 
are significant differences in the three flows at this 
position, and also the statistical properties have been 
carefully documented. 

Before progressing further, a new term must be intro- 
duced and it is the line of maximum maximorum of the 
correlations . Physically what is meant is that a line 
is drawn through the point on each correlation level 
which is separated the farthest amount from the reference 

point (i.e. (5 2 + ? 2 )^ or (5 2 + ? ) 2 is a maximum) . 
12 2 3 

With this line, it is then possible to get an idea of the 
flow pattern. 

In Figure 50, the isocorrelation curves in the x-z 

plane for the case of the unconfined, coflowing jet are 

shown. The curves are very elongated in the streamwise 

direction and close to being symmetric with respect to 

the ? axis. This is not the case about the £ axis, 

1 3 

which is in contrast with contours presented for a fully 
developed turbulent boundary layer. 1 } The difference 
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presumbly is due to the fact that part of the flow field 
lies in the potential core while the rest can be consider- 
ed a mixing region. The line S, marking the points of 
maximum maximorum of the correlations, is approximately 
straight as expected. 

The isocorrelation curves in the x-z plane for the 
case of' the turbulent jet flowing over the plate are 
shown in Figure 51. Once again, the curves are elongated 
in the streamwise direction. There quite clearly is 
not any axis of symmetry. The pattern formed by the iso- 
correlation curves is very complex and is composed partly 
of two different lobes of high correlation value levels. 

The line S is drawn through the points of maximum maximorum 
of the correlations as was done in Figure 52. Downstream 
of the fixed reference point, the line is convex towards 
the surface of the wall. Upstream of the reference 
point, it is slightly more difficult to construct the 
line but it also seems to be convex towards the plate's 
surface. Sternburg^^ has offered a possible explanation 
of a qualitative nature for the observed shape of the 
S-lines. He theorizes that this may be caused by the 
mean shear, deforming the turbulence flow pattern as 
determined by its vorticity distribution. The maximum 
correlation will be determined by a point of the vorticity 
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pattern that passes the fixed laser control volume at t 
equal to zero and is observed at some time, x, later. 

The larger eddies which have a longer correlation distance 
will rotate due to the distortion of the mean shear flow 
in a clockwise sense and with a center of rotation in 
x-z plane. As the second, movable control volume ap- 
proaches the plate, a point of the flow pattern passing 
the fixed laser "probe" at x=0 will be displaced on the 
average from the wall by this rotation for x > 0, while 
for x < 0 this point has arrived from an upstream position 
also at a greater distance from the wall. The rough 
picture would lead to an S-line convex towards the wall 
below the center of rotation, a fairly straight line near 
the center and it will bu curved concave to the wall 
above the center. 

The isocorrelation contours in the x-z plane ^f or 

the flow of the turbulent jet over the flap are presented 

in Figure 52. As is the case for the plate configuration, 

the pattern formed by the isocorrelation contours is clearly 

not symmetric about either the £ or the ? axis. Before 

1 3 

it has been observed that the previous two iso-contours de- 
monstrated an elongation in the streamwise direction. 

For the flow over the flap, the elongation exists but the 
pattern formed is rotated- a sizeable 'amount from the 
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longitudinal direction. Through the points of maximum 
maximorum of the correlations/ the line S is drawn. Before 
applying Sterburg’s approach/ -it is first necessary to re- 
call the geometry of the flap and the relation position 
of the reference "probe”. In Figure 7(a)/ a view of the 
jet/flap assembly is shown. After a distance of approximate- 
ly four diameters downstream from- the exit plane of the 
nozzle, the tangent to the surface of the flap rotates quick- 
ly away from the horizontal direction. Thus, any position 

for the second control, volume downstream of the x/2r = 4 

o 

plane would be over the curved section of the flap ' s 
surface. Now, consider the shape of line S in Figure 47. 
Downstream of the reference point, the line heads down 
towards the wall surface -.and is also convex to the flap. 
Upstream the line is also inclined with a negative slope 
from the x-direction, and is essentially straight. Clear- 
ly, it is possible to say from .observing the iso- 
correlation contours that the flow is being turned by the 
flap. Secondly, considering the change in shape of the 
flap surface and the convex shape of the line S down- 
stream, Sternburg's ideas can be reasonably applied once 
again. For the clockwise rotation of the large eddies, 
a point of the flow pattern passing the reference control 
volume will tend to both be turned downward by the flap ’ s 
surface and also be displaced on the average from the wall 
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by this rotation for a positive delay time. The flow 
pattern upstream of the fixed reference point indicates 
the direction from which the flow is coming. 

In Figure 53, the isocorrelation curves in the y-z 
plane for the unconfined coflowing turbulent jet flow field 
are shown. The eccentricity of each of the curves does 
vary but is essentially zero and thus the isocorrelation 
levels are close to being circles in shape. Hence, 
symmetry exists about both the and the axes. Though 
Sternburg's ideas cannot be applied in the y-z plane, the 
line S can be constructed and from its shape, a rough 
estimate- of the flow pattern can be obtained. For the 
unconfined jet, the line of maximum maximorum can be 
drawn, starting from the center of field, off in any radial 
direction. One line is arbitrarily constructed for the 
sake of clarity. From observing the line S and the iso- 
curves, an indication that the jet's flow field is spread- 
ing out equally in all directions is apparent. The lower 
correlation, levels farther out from the center line re- 
present fluid particles which have been separated from 
the potential core the longest period of time. 

Consider next the isocorrelation curves in the y-z 
plane for the flow of the turbulent jet over the plate 
Figure 54. The change in the shape of the curves caused 
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by tiie presence of tlie plate is apparent. The iso- 
correlation curves are elongated in the lateral direction. 
Although the flow pattern is clearly not symmetric about 

the £ axis, symmetry does exist with respect to the £ 

* 2 3 

axis. As was the case in the x-z plane, the line S is 

again somewhat more difficult to construct. Allowing 

some author's discretion, the line S moves out from the 

center of the flow in a slightly convex manner with 

respect to the plate's surface. Two impressions can be 

drawn. Initially, the- turbulent jet is drawn downwards 

towards the plate. Secondly, there exists an indication 

that the flow -may be rolling up in the outer regions, 

or at least some type of canted rotating motion of the 

larger eddies to the longitudinal downstream is being 

observed. It is instructive to compare the differing 

flow patterns displayed in Figures 54 and 55. The width 

of area enclosed by the isocorrelation curves is 

significantly smaller for the flow over the plate for 

£ >0, about the same for £ =0 and somewhat larger 

3 3 

for the £ <0. The second observation is -reinforced by 

3 

the measurements made of the mixing widths and the 
intermittency profiles for the two flows. 

Lastly, the isocorrelation curves in the y-z plane 
for the flow of the turbulent jet over the flap are shown 
in Figure 55. The' isocorrelation curves are markedly 
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elongated in the z direction. The flow pattern is 

symmetric about the ? axis but asymmetric with respect 

* 3 

to the 5 axis. The- line S is drawn appropriately 
2 

through the flow field pattern. The first impression is 
one of noticing the narrowness of the flow pattern. As 
has been mentioned previously, the flap is much more 
effective in amplifying the turbulence levels of the tur- 
bulent jet. This more rapid and intense mixing action 
would manifest itself in a decrease in the coherence of 
the flow field pattern. Thus, in the outer mixing region, 
the turbulent jet is more apt to have forgotten its 
previous flow condi tons, and correlations between the 
center of field and a second point farther out in the 
mixing layer would be weaker. Secondly, consider the 
direction of the line S . Downstream of the fixed re- 
ference point, the line is clearly heading in a negative 

5 (i.e. - z) direction. Upstream,- there is some 

3 

arbitrariness involved but the line also has a negative 

£ heading. From this flow pattern the impression is 
3 

made again that flap is quite effective in turning the 
flow. At a positive delay time, the particles that were 
at the reference point have moved steadily downwards to- 
wards the surface of the flap, and thus a stronger cor- 
relation level exists. The same type of agreement can 
also be made for the flow .upstream of the reference point. 
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Section 4 . 8 Convection Velocities 


(c C \ 

The possibility exists of defining many kinds 
of convection velocities. Favre, Gaviglio and Dumas 
showed that convection velocities can be determined from 
iso- space- time velocity correlation curves R , (£,0,0, T) . 
Favre et. al. then defined the convection velocity as 
follows : 


Ucon = £ /t 

1 max 


where t is the delay time for which R , has a maximum 

iud - ■“ 


x 1 


value . 

In particular in the case of relatively strong in- 
homogeniety taking in the dimensionless time U t/L, an 
average value of U across the distance £ , i.e. 

1 ,5- , 

0=7- J rt i d£ 1 U (x + £•;) = [U (x ) + U (x + £ ) ] 

1 i° 111 1 z 11 l 1 


where L is any suitable length scale. 

The convection velocities determined for the three 
different flow configurations are shown in Table 2. Here 
for a given longitudinal displacement, nondimensionalized 
by twice the radius of the jet, between the two measuring 
laser "probes", the non-dimensionalized delay time and 
the ratio of the convection velocity to the mean velocity 
are given. Mean velocity here refers to the average 
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TABLE 2 


Convection Velocities in Various Flows 


Coflowing Jet Coflowing Jet/Plate 


5 1 

U 

1 T 

Ucon 

£ 

1 

U . 

1 <r 

Ucon 

2r 

o 

2r 

o 

°i 

2r 

o 

2r 

o 


- 1 752 

i 

-.561 

' ■ 1.34 

-.761 

-.786 

.968 

-.595 

-.527 

1.12 

-.595 

-.587 

1.01 

0 

0 

1 

0 

0 

1 

.177 

.210 

.843 

.262 

.262 

1 

.610 

. 606 

1.01 

.476 

.393 

1.21 

.942 

.921 

1.02 

1.14 

1.05 

1.08 

1.22 

1.17 

1.04 

1.42 

1.28 

1.11 

1.59 

1.51 

1.05 

2.23 

1.79 

1.25 

1.82 - 

1.65 

1.10 

3.06 

2.42 

1.26 

2.73 

2.49 

1.09 

5.09 

3.07 

1.56 


S 

l 

u- 

1 T 

Ucon 

2r 

o 

2r o 


-1.49 

1.28 

1.16 

-1.03 

-.925 ' 

1.11 

-.480 

'-.40 

1.20 

0 

0 

1 

.426 

1.76 

.242 

.770 

2.98 

.258 

1.14 

3.83 

.297 

1.86 

6.26 

.297 

2.55 

8.50 

.300 
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value* determined using Favre's definition. 

Focussing initially on the convection velocity 
calculated for the unconfined, coflowing turbulent jet, 
the ratio of the Ucon to IT is very close to unity, and, 
indeed, remains fairly constant. Thus, in this case,' at 
least, the convection velocity appears to be almost in- 
dependent .of the separation distance. Next, consider the 
case of the turbulent jet flowing over the plate. Here 
someth ng quite interesting happens. As the flow pro- 
gresses downstream, the convection velocity has ap- 
proximately the same value as the mean velocity for 
?^/2r Q < 0 but for £^/2r Q increasing positively, the 

ratio Ucon/U is greater than unity. This phenomenon also 

^ ( 6 5 ) 

occurs in turbulent boundary layers over flat plates. 

To explain the reason for this behavior, recall the 

vertical height above the surface of the plate at which 

these measurements are taken and also the shape of the 

mean velocity profile for this flow in the center vertical 

plane. Notice that there is a region of faster moving 

fluid particles closer to the wall surface. Also, as 

the flow progresses downstream, the profile approaches 

that found for a classical wall jet. Hence though the 

ratio of Ucon/U^ may be greater than unity, the ratio of 

Uc/U max will never be. 
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Finally, consider the flow of the turbulent jet 
over the flap. Notice that for £ /2r Q < 0, the value 
of the ratio of Ucon/U^ is slightly greater than unity. 
This can be explained using the same logic that is used 
in the previous paragraph. The reason that Ucon/U > 1 
occurs sooner than is the case for the flow over the 
plate is that the mean velocity profile over the curved 
surface in the center vertical plane is transformed much 
more effectively and rapidly from the ideal jet's "flat 
top" to the wall jet shape. For 5^/2r Q > 0, the con- 
vection velocity has decreased substantially. Pre- 
sumably this is due to the fact that the flow is under- 
going a change in direction. 

Section 4 . 9 Pressure-Velocity Correlations 

Additional information concerning the turbulence 
structure of the various flow fields can be gained from 
measurements of the pressure fluctuations at both the 
wall and in the turbulent jet and correlating those 
signals with the fluctuating turbulent velocities in 
the potential core and in the shearing region. . 

With the orgin of the coordinate system at the 
location of either the surface port or the pressure 
probe, the follwoing space-time correlations are mea- 
sured: 
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P (O/O/t) U (£, /? , *t + t) ' 

2 3 

^ / (x,y, z,t, £ ^ £ / 5 / * t ) ~ 

|pj ^ 2 g —* "‘'‘ 1 ■ ■ n 

[p(o,o;t) 2 K [u(? ,5 jt+t) 2 ]" 2 

2 3 


where p is the static pressure measured at the wall or 
in the flow field. 

Since the primary focus of this segment of the in- 
vestigation is to determine the relationship between the 
pressure and velocity fields, the coherence is plotted 
for various pressure and velocity monitoring locations. 
Recall that coherence is between two signals defined as 
follows : 



12 G G 

11 22 

where G , G are the Fourier Transforms of the in- 
11 22 

dividual autocorrelation functions and G is the Fourier 

12 

Transform of the cross correlation function. 

Case I Wall Pressure Fluctuations 
Cross correlations between fluctuating pressure 
signals measured at surface ports in the plate and flap 
and turbulent velocity signals monitored at various 
locations in the flow field are determined. For the flat 
plate case, coherences are determined with the laser 
probe at 5 /2r = 0.2 with £ /2r and £ /2r set equal 
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4 are 


to zero. Examples of the results for x/2r Q = 
shown in Figure 56. Here the lateral positioning of the 
pressure port and laser control volume varies from |y[/2r 
0 to |y'f/2r 1. It is found that inside the lip of the 

nozzle, the velocity spectra and pressure spectra are 
characterized by large peaks at, in this instance, ap- 
proximately 3G0HZ. The coherence between the pressure 
and velocity signals is clearly the strongest at this 
peak frequency. As the pressure and velocity monitoring 
"probes" are moved outward from |yj/2r Q =0, the spectral 
peaks are attenuated and the resultant decrease in the 
coherence is noted in Figure 56. 

Coherences are also determined for the flow over 

the plate with £ /2r varying from 0 to 1, £ /2r 

2 o l o 

held equal to zero and £ /2r set equal to 0.5 (Figure 

3 ° 

57) . Physically, the laser control volume is held 
fixed at the location x/2r Q = 4, jyJ/2r Q = 0, and z/2r Q = 
0.5 with the fluctuating velocity signal being correlated 
with surface pressures measured at x/2r Q = 4, z/2r Q = 0 
and |yj/2r = 0,0-5 and 1. It is found that the velocity 
spectra is^ not as peaked as is the case at z/2r Q = 0.2 
and, hence, the coherence between the pressure and. 
velocity signals is diminished. Also, the coherence 
decreases as the surface probe is moved laterally out- 
ward from the jet centerline. 
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For the flow over the flap, coherences are deter- 
mined between the pressures monitored at the surface and 
the velocities at the |y[/2r Q = 0 location. ’ The downstream 
location is held fixed at x/2r Q = 4. In Figure 58, 
coherences are shown with the laser "probe" held fixed 
at |y|/2r Q = 0 and z/2r Q = 0.2. while in . Figure • 59 , the 
laser measuring volume is located at |y|/2r Q = 0 and 
z/2r Q = 0.5. Many of the' same observations that were 
made concerning the flow over the plate can be made 
once again. The velocity spectra is found to be much 
more peaked closer to the surface of plate. The coherence 
is larger in magnitude between the wall surface pressure 

and the velocity fluctuations with £ /2r = 0.2. Also, 

3 ° 

the coherence is the strongest at approximately 300Hz. 

Case II Pressure Fluctuations at Probe 

Cross correlations and coherences are determined- between 
the static pressure measured by a static pressure probe 
and the turbulent velocity fluctuations measured at the 
jet axis. The pressure is monitored at three lateral 
locations, (i.e. j;yj/2r Q = 0, 0.5, 1.0) with the velocity 
control volume remaining fixed. The vertical location 
for the pressure probe and LDV volume remains fixed at 
z/2r Q = 0.5. 

Figure 60 shows the change in the value of the 
correlation coefficient at a given lateral location as 
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Figure 58 Wall Pressure and Velocity Coherence. 5 /2r = 0, 

/O 2 O 

? 3 // r o = '02. (Xj = 5.1, DB reference to Dynamic 
Pressure) 


147 


Figure 59 Wall Pressure and Velocity Coherence. £ /2r = 

£ /2r ® 0.5 (A = 5.1, DB referenced to dynamic 

O'-' J 

Pressure) 





the flow progresses downstream in each of the three con- 
figurations. The decrease in the correlation coefficient 
as x/2r Q increases is. apparent. Also, the flow of the 
unconfined, coflowing turbulent jet maintains the highest 
value of the coefficient at each measurement position. 

For the unconfined, coflowing jet, the coherences be- 
tween the fluctuating pressure and velocity signals are 
shown in Figure 61 for the downstream location, x/2r Q = 4. 
It is found that both the pressure and velocity spectra 
are characterized by large peaks at approximately 300Hz . 
This clearly is the frequency for which the coherence 
is- a maximum. As is evidenced in the case of wall pressure 
fluctuations, the coherence between the pressure and 
velocity signals decreases as the pressure probe is 
moved laterally outward. 

For purposes of comparison, the pressure and velocity 
spectra, and coherence at x/2r Q = 8 at the jet axis. are 
also shown in Figure -62. The peak which characterizes 
both the pressure and velocity spectra is absent and the 
coherence has been reduced significantly. At this dis- 
tance downstream from the exit plane, the turbulence has 
become much more fully developed. 

Next, consider the case of the turbulent jet flow- 
ing over the flat plate for the downstream location 
x/2r Q = 4 (Figure 63) . The pressure spectra is again 
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Figure 63 Static Pressure and Velocity Coherence. 5^/2r Q = 0, 

? /2r Q = o. (*• = 5.1, DB referenced to Dynamic Pressure) 

3-1 






found to be markedly peaked at 300Hz at |yj/2r =0 and 

o 

this "bump" is both broadened and attenuated out from the 
jet axis. The velocity spectra though not peaked' as it 
is in the case of the unconfined jet is not at all similar 
to the fully turbulent example at x/2r Q = 8. The coherence 
of the two fluctuating signals has been significantly re- 
duced with the plate in the flow field, and it decreases 
with increasing lateral separation of the pressure probe 
and laser control volume. The coherence obtained from 
the flow field of the turbulent jet and flap at the 
downstream location x/2r Q = 4 are shown in Figure 64. 

The peak in the velocity spectra has been virtually 
eliminated. . It is interesting to compare the coherence 

of each £ /2r position for this flow with the jet/plate 

■ 2 o 

and unconfined jet flow fields. The coherence for the 
two confined jet flows are quite close in magnitude. The 
transformation of the turbulent jet into- either a classical 
wall jet or curved wall jet does decrease the relationship 
between the fluctuating pressure and velocity fields . 



Chapter 5 


SUMMARY 

A summary of the significant results detailed in 
this experiment will now be given. The presence of the 
confining surfaces was found to be of importance in 
the resultant effects on both the mean and fluctuating 
part of the turbulent velocity field of the jet. 

(1) The mean velocity was measured for each of the 
three flow configurations at the jet axis (z/2r Q = 0.5) 
at different downstream locations . The decay of the mean 
velocity of the turbulent jet was accelerated by the 
presence of the confining surfaces. The, flap was also 
found to be the most effective in diminishing the x- 
directed momentum. For an increasing value of X ^ , the 
effectiveness of the flap in diminishing the x-directed 

•momentum decreased. Evidence existed that as X ^ ap- 
proached infinity, the flow would not remain attached. 

(2) The flap seemed to break up the top hat 
characteristic shape of the mean velocity profile at 
the exit plane much more rapidly. The width of the 
flow field increased at a faster rate for the flow over 
the flap. The parameter X^ also influenced the width of 
the mean velocity profiles . For the case of the turbulent 
jet/flap, an increase in Xj caused a resultant decrease 
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in the mixing width, y . 

(3) It was found that for the flow over the plate, 
the mean velocity field tended toward the fully develop- 
ed turbulent shape in different horizontal planes faster 
than in the vertical plane at y/2r Q = 0. Also, the flap 
seemed more effective in changing the flat velocity profile 

characteristic of a circular nozzle into a wall jet profile. 

(4) The confining surfaces caused a more rapid in- 
crease in the turbulent intensities in the x-direction 
at the jet axis. Here, as well, the value of A ^ was 
also found to be important. At the large value of A^, 
the turbulent intensities for the flow over the flap 
increased at a significantly less rapid pace than was 

the case for A. equal to 5.1. An indication that the 
J 

flap and the plate serve to break up the potential core 
more rapidly was also noted and discussed. 

(5) Turbulent intensities in the lateral, y, 
direction for the three flow configurations were also 
measured. The turbulent intensity was found to be much 
greater in the shearing layer of the unconfined jet than 
in either of the other two configurations . Some evidence 
thus exists which . indicate a damping of the turbulence 

in the y direction by the presence of the confining 
surfaces . 
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(6) Autocorrelations of the longitudinal velocity 
fluctuations were taken at the jet axis and at the lateral 
position corresponding to the lip of the nozzle at various 
x/2r Q locations. .The autocorrelations taken in the potential 
core regions were characterized by a damped sinusoidal 
nature. The dominant frequency in the autocorrelation curve 
is theorized to be due to the shedding of vortices (i.e. 
large sca'le structures) at the exit plane of the nozzle. 

As the vortices pass downstream, the associated fluctuating 
pressure field causes velocity fluctuations which are 
monitored. The effects of the plate and especially the flap 
on the , autocorrelation were apparent, causing a shift and 
attenuation of the frequency at z/2r Q = 4. The changes in 
the autocorrelation curves as the flow progresses down- 
stream are also noted. In the mixing layers of the flow 
fields, the dominant frequency shape was absent due to the 
large increase in the turbulent level and it was very difficult 
to determine any effects caused by presence of the confining 
surfaces . 

(7) Longitudinal integral length scales and Taylor 
microscales were determined. The integral scales for the 
flap were seen to have a rapid increase and then decrease 
in magnitude as the flow approached the start of the curved 
portion of the surface. For the unconfined turbulent jet, 
the growth of the integral scale at |y|/2r Q was linear, as 


158 



predicted. At best, the integral scale development for the 
three different configurations was a weak function of Aj. 

The value of A ^ did strongly influence the magnitude of the 
Taylor microscale, A^, at the lateral position corresponding 
to the lip of the nozzle. 

(8) Power spectra for the different confining surfaces 
reflected the change observed in the autocorrelation curves. 

At x/2r Q = 4, the flows of the unconfined jet and jet/plate 
were still peaked in nature at the jet axis but for. the flow 
over the flap, this "bump" was absent. The "bump" or peak 
in the power spectra corresponds to the passing downstream 
of the large scale vortex rings. As the flow becomes more 
fully developed due either to the amplification of the 
turbulence by the confining walls or simply the flow pro- 
gressing downstream, the relatively small velocity fluctuations 
due to the shedding vortices become less apparent and 
essentially become "lost" in the increased turbulent levels. 

It was interesting to note that with regard to the flap 
configuration, the absence of the peak was true for A ^ = 

5.1 and A ^ = 10.88. This is inspite of the fact that for 
A j = 10.88, the flow had not decelerated as much nor had 
there been as large an increase in the turbulent intensities 
in the central part of the velocity field. The downstream 
development of the power spectral densities at' the jet axis, 
[y|/2r Q = 0, and out at the lip of the nozzle, jyJ/2r Q = 0.5, 
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were also presented. At the lateral position downstream 
from the lip of the nozzle, the power spectra consistently 
increased the low frequency content of the spectra as the 
flow progressed downstream. 

(9) Auto and cross correlations were determined for 

the concentration field at various downstream locations. 

At x/2r = 4, the autocorrelation coefficient curve was 
o 

found to possess a damped sinusoidal nature for the three 

different flow configurations. The concentration fluctuations 

at y = ± 0 . 5 were found to be related both at x/2r = 4 

o 

and x/2r Q = 8. The maximum absolute values of the cross 
correlation coefficient did decrease sizably from four to 
eight diameters downstream of the exit plane, however. The 
power spectra for the jet, jet/plate, and jet flap all ex- 
hibited a peak at x/2r Q = 4 at approximately the same • frequency 
as that found in the core velocity spectra. Here again, 
evidence exists in reinforcement of the vortex model of the 
jet. The plate and the flap did seem to diminish the peak 
and also seemed to cause a shift of energy to lower fre- 
quencies in the spectra. 

(10) Space-time correlations for the turbulent velocity 
fields were obtained, and the lines of maximum maximorum 
were constructed for both the x-z plane and the y-z plane. 

For the unconfined jet, the iso-correlation levels were 
symmetric with respect to the axis though very elongated. 
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The line of maximum maximorum was straight and parallel to 
the £ axis. In the y-z plane, the isocorrelation levels 
were circular in shape about the E, ^ axis. For the flow over 
the plate , the isocorrelation curves in the x-z plane were 
very complicated and elongated in the £ direction and the 
line of maximum maximorum was curved convex towards the 
wall both upstream and downstream of the fixed reference 
"probe". For the flow over tlie flap the isocorrelation 
curves were elongated and rotated down towards the surface 
in the x-z plane, evidence also existed as to confirm the 
effectiveness of the flap in turning the flow. The iso- 
correlation contours for the turbulent flow fields 
demonstrate the existence of large scale structures. Con- 
sider the fact that the shape of the contours is significantly 
different for each of the three flow configurations in both 
the x-z and the y-z planes. The flow "knows" whether or 
not a confining surface is present and also "senses" if the 
wall is curved or flat. This "knowledge" is then seemingly 
transmitted throughout the flow field. 

(11) Convection velocities were determined for the 
three flow configurations. The convection velocity deter- 
mined for the unconfined jet was found to be approximately 
equal to the local mean velocity and also to be independent 
of the separation distance between the "probes". For the 
cases of the flows over the plate and flap, the convection 
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velocity was greater than the local mean velocity in several 
instances when the local velocity was considerably less than 
the maximum velocity at that particular downstream location. 

(12) Pressure velocity correlations using both the 
static pressure probe and the surface ports yielded strong 
evidence that as the flow progresses^ downstream, and the 
flow becomes a fully developed turbulent flow, the relation- 
ship between the pressure and velocity fields diminishes. 

For the first several diameters downstream from the exit 
plane when the pressure and velocity spectra peak at 
approximately 300 Hz, the coherence between the two 
fluctuating fields is the strongest. 

This experimental investigation has shown that large 
scale structures do exist in the near field of the turbulent 
jet. The two point velocity and concentration measurements 
document this observation. The plate and the flap were 
found to strongly influence the large scale structure of 
the jet. 

Several conclusions can be reached. First, the flap 
serves to diminish the extent of influence of the large 
scale structures in the flow field. This should help in 
reducing the structural loads the flap or wing must endure. 
Instead of one dominant frequency in the velocity field, the 
turbulent fluctuations in the entire power spectrum would 
be amplified. The implication would then be to mount the 
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nozzle as close as possible to the curved portion of the 
flap. The second conclusion is in reference to comparing 
the velocity flow field with the two different velocity 
ratios, X ^ . There exists considerable doubt as to whether 
or not the jet flow will even remain attached as the 
velocity ratio between the inner jet flow and the tunnel 
flow tends towards infinity. The significance of 
measurements made in a static testing configuration 
( X j -> <°) is quite questionable. 
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APPENDIX A 


LDV Parameters and Velocity Field Properties 

The LDV parameters include such items as the 
characteristics of the data processing unit, the control 
volume, size and the signal-to-noise ratio of the Doppler 
bursts. Initially, consider the effects of the probe 
volume size. For the optical arrangement used in this 
investigation, the control volume was estimated to be the 
order of 0.05 mm (diameter) x 1 mm (length) . 

Mean Velocity Analysis for Turbulent Flow 

The mean velocity, Uam, determined by a finite sized 
probe includes a space as well as a time averaging. Con- 
sider the following expression: 

gam . JL //[(l/T)/Jj(A,t)dt]dA (1) 

where U(A,t) is the instantaneous streamwise velocity, 
t is time and T is time interval for measurement. 

For simplicity, assume a rectangular probe volume 
and a two dimensional flow field then 

Uam = (l/w)/^[(l/T)/^U(y,t)dt3dy (2) 

where w is probe volume dimension in the direction of ' 
the mean, velocity gradient, y, and y equals y-y Q . 
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Next, let U(A,t) vary linearly with y across the 
probe volume or 

U (A, t) - U (t) + Sy 

Then U(A,t) = U (A) = U 4- Sy (3) 

c 

Substituting back into (2) 

Uam = (1/w) 

Thus Uam = U (4) 

c 

The, conclusion being that the summation over a 

finite volume and time yields a go*od estimate of the mean 

velocity at the center of the probe when the average 

fluid velocity is linear. 

Root Mean Square Velocity Measurement 

The spatial and temporally averaged mean-square 
2 

fluctuation, (u^) of the streamwise fluid velocity 
flowing through the probe volume is given by 

(U m )2 = (1/A)//E(1 /T)/.q (U (A, t) - U c ) 2 dt]dA (5) 

when the velocity fluctuations are referenced to the 
time average velocity at the center of the probe volume. 

If we again assume that U(A,t) *» U{y,t), a function 
of y alone and also that the volume is rectangular then 

<V 2 = (1/w) lw/2 [( 1 /T)/o (u (y,t) - U c ) 2 dt]dy (6) 
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Let us next assume that the velocity fluctuations 
are inversely proportional to the mean velocity. This 
makes sense, as an approximation in the mixing region of 


the- jet. 

Then u* Cy) ' = K/U(y) = K/{5 + Sy) (7) 

v 

Therefore 

u 1 (y) 2 = (t/T)S^ [U (v, t) - U(y)] 2 dt (8) 
u* (y) 2 » K 2 /(0 + Sy) 2 (9) 

Rearranging ( 6 ) 

u- m 2 = (1/w )f ™/2 E(Vtt)/-; “ u(y)) 2 dt 

+ (O(y) - U c > 2 ]dy . (10} 

Using (7) and (9 ) , (10) may be integrated to yield 

2 = S 2 w 2 + 4Kf (11) 

^ 12 4U c 2 - S 2 w 2 -• 

But K 2 /U 2 is the mean square fluctuation at y = 0 
c 

and therefore 

• (us ) 2 ^ §l w . 2 .. + 4 K 2 : (uq :_) 2 (12) 

m 12 4K 2 - S 2 w 2 (U Q . * ) 2 

The first term on the right hand side of (12) re- 
presents the "false turbulence" seen by the LDV due to 
a mean velocity gradient across the control volume. The 
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second term is a combination of the actual turbulent 
velocity variance and "false turbulence" due to change 
in mean square fluctuations across the probe volume. 

For the experimental investigation which has been 
reported here, let us make an estimate of the magnitude 
of the "false turbulence" measured by the LDV, under the 
given assumptions. 

Typical values of the parameters would be the 
following: 

u , 

~ 0.1 w = 1 mm K — 1 x 10 2 m 2 /sec 2 

o 

= 10 m/sec S = 1 x 10 3 m/sec/m. 

Then, the actual variance of the flow field is: 

(u Q ' ) 2 “1 m 2 /sec 2 

The measured velocity variance is: 

(u_ ' ) 2 = 
m • 


( u ra ) 2 ' = 
= 


[ ( 1 x 10 3 ) 2 (1 x 10~ 3 ) 2 + 

12 

4 W.) tl.) , m 2 /sec 2 

4 (10 4 ) - (3 x *10 3 ) (1 x .10 3 ) ] 

.083 + 1.00 m 2 /sec 2 
1.083 m 2 /sec 2 
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From the very simplified analysis , there are several 
conclusions that can be made. First, both the size of 
the control volume and the magnitude of the mean velocity 
gradient can have a non-negligible effect on the LDV 
turbulent intensity measurements. If at all possible, 
the smaller dimension of the control volume (length or 
height) should be aligned parallel to the direction of 
the largest mean velocity gradient. 

Secondly, for the assumed form of u' (y) , the velocity 
variance caused by the change in the mean square fluctua- 
tions across the probe volume is negligible. If a linear 
relationship between the mean velocity and fluctuating 
velocity were proposed (i.e. u' (.y) = KUfy)) taken the 
contribution to the total variance would take the 
following form: 


S 2 w 2 A 2 

12 

where S, w would take on the same value as before but 
K = 1 x 1(T 2 . 

Different functional forms for u‘ (y) where assumed 
and the results are shown in Table 3. 

Next, the characteristics of the signal processing 
unit (in this case a phase-locked loop) are also of 
paramount performance in the determining the accuracy of 
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TABLE 3 


Evaluation of Effect of Control Volume 
Size on Mean Square Fluctuations 
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the LDV measurements. The phase locked-loop used is able 
to respond at a 14kHz rate to deviations from the 
centered or mean frequency (velocity) up to ± 28%. 

The unprocessed Doppler signal .was fed into a Bio- 
mation Transient Recorder Model 8100. The Transient 
Recorder was used as it allows different segments of the 
Doppler signal to be captured and displayed on an . 
oscilloscope. From this visual display, the frequency 
of the individual Doppler bursts can be determined, and, 
hence, if enough samples are displayed, a histogram 
can be constructed. Two histograms were obtained, one 
in the potentional core region of the freely expanding 
jet flow and the other determined in the shearing region 

of the velocity field; For the flow in the potential 

\ 

core region, the ± 28% range was more than adequate to 
handle the small turbulent fluctuations. In the shearing 
regions, though the turbulence was quite high (i.e. 

(u 2 ) /U o > 0.1), the phase locked loop's range still 
proved to be quite adequate. Out of the approximately 
500 instantaneous frequencies measured, only about 
5% exceeded the limit of ± 28% of centered frequency. 

The signal-to-noise ratios for the Doppler signal 

\ 

were determined in a variety of ways. The first signal 
to noise ratios determined were for the signal before 
it was fed into the phase locked-loop. The filtered 
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and amplified Doppler signal was analyzed by use of a 
Nelson Ross Spectrum Analyzer. The analyzer displayed 
the frequency content of the incoming signal. Here, two 
approaches were taken. First, with the Doppler signal 
present, an estimate of the signal to noise ratio could 
be determined from the visual display of the analyzer 
an' oscilloscope. This ratio varied from about 100 to 1, 
for the photodiode in a dual beam forward scattering 
mode, 75 to 1 for the reference beam mode (diode) , to 
approximately 20 to 1 for the photomultiplier tube in 
a back scattering location. Another estimate was made 
with the level of the Doppler signal compared to the 
noise present when one of the crossing laser beams is 
blocked. For the photomultiplier tube, the d.c. current 
was increased so as to be comparable to that present 
when the beams intersect in the- test section. The signal 
to noise ratios determined from this procedure varied 
from about 100 to 1 for the dual beam mode (diode) , 75 
to 1 for the reference beam mode (diode) , to about 15 
to 1 for the photomultiplier tube in a back scattering 
mode. Since only the zero crossings are of interest, 
it would seem that the presence of the Doppler signal 
seems to override the noise in the case of the photo- 
multiplier tube. The noise is thus superimposed on the 
Doppler signal. 
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The third approach taken to determine the signal to 
noise ratios for the LDV setup involved using the output 
voltage, E, of the phase locked-loop. Here the Erms of 
the Doppler signal was compared to the Erms output of the 
phase locked-loop - when, again, one of the crossing laser 
beams was blocked. The ratios determined in this manner 
ranged from 50 to 1 for the forward scattering dual beam 
mode (diode) , 20 tol for reference beam mode . (diode) , to 
about 15 to 1 for the back scattering mode using the 
photomultiplier tube. 

With the wide variance of signal-to-noise ratios 
for the forward scattering modes as compared to the 
back scattering mode, a logical step would be to compare 
the data obtained from the two setups. A sample of the 
results of this comparison is shown in Table 4. 

In the sample data presented in Table 4 as well as 
for the remaining comparisons made the mean frequencies 
measured by all three modes of operation are extremely 
close. The rms voltages (i.e. turbulent intensity, 
(u 2 5%/Uo) are nearly identical for the two forward 
scattering modes. Erms obtained from the photomultiplier 
tube used in a back scattering configuration is slightly 
less, though still reasonably close. This trend was 
consistent for all the comparisons that were made, both 
in the relatively quiet potential core part of the flow 
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TABLE 4 


Case 1 


(Potential 

Core) 


Case 2 


(Shearing 

Region) 


Mode 

Mean Frequency (MHz) 

Erms (Volts) 

Dual 

2.246 

.460 

Reference 

2.243 

.465 

Back Scatter 

2.247 

.435 


Mode 

Mean Frequency (MHz ) 

Erms (Volts) 

Dual 

1.356 

1.08. 

Reference 

1.358 

1.05 

Back Scatter 

1.343 

0.99 
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as well as in the highly turbulent shearing region. 
Presumably, this slight reduction in the turbulence 
perceived by the LDV in the back scattering mode is due 
to the fact that light scattered in this direction is 
several orders of magnitude less than is the case for a 
forward scattering mode. Even with the higher potential 
gain possible the photomultiplier tube is, with the 
added constraint of a low noise level, not capable of 
detecting 1 all the particles that are passing through the 
control volume. In essence, the contribution of these 
particles is "lost" in the noise. 

Along with. this discussion of signal-to-noise ratios 
it may prove to be instructive to compare in a qualitative 
manner, as well, the Doppler signals observed by the 
three different modes of LDV operation . (Figure 64). These 
drawings were made by use of the Biomation Transient 
Recorder. The control volume was located in the fully 
developed turbulent region of the freely expanding co- 
flowing jet. There are two distinctions that can be 
made. One, the Doppler "bursts" in the forward scattering 
modes are fairly uniform envelopes. With the back 
scattering mode, however,, there is quite a bit more 
ampltitude modulation of the signal. The second dis- 
tinction deals with a comparison of the "dead" time 
(i.e. no particles in the scattering volume) . For the 
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Volts 



Figure 65(a) 


Doppler Signa 





forward scattering modes, the "dead" time is distinct 
and easily distinguishable for the Doppler bursts. Such 
is not the case for the back scattering mode. It appears 
that some frequency component is virtually present the 
entire interval shown. The sharp distinction between 
signal and dead time is not longer present. 
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APPENDIX B 


Particle Dynamics 

The LDV measures the particle velocities and the 
flow velocity is then inferred from this information. 

The question of how closely the particles follow the 
actual flow field is an extremely important one when 
considering the accuracy and reliability of the LDV. 

Consider the equation of motion of a particle in a 
non-uniform flow under the assumption that the only 
effective force acting on the particle is the aerodynamic 
drag, and that the gas is much less dense than the 
particles (pg/pp << 1) then: 

— • P = 3/4 * ^ Re .* C • (U - U ) (51) (1) 

dt pp D 2 P D * ? 

where Ug 
Up 

Re P 


= velocity of gas 
= velocity of particle 



y 


and C_ = f (Re , Mach no., Knudsen no.) 
D p 

Using Stokes Law: 



( 2 ) 
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Inserting (2) in (1) 


^U£ = 3/4 -2i — Rep (D - U ) (3) 

dt Re p D 2 g P 


where K, 


^| = K (Ug - Up) 

the Stokes coefficient is given by K = 


(4) 


18M 


p p d ‘ 


For this investigation, 

D = . 3ym - .3 x 10*" 6 m 


yair = 1.796 x 10" 5 kg/m sec 

eg = 1.4 or pp = 1.4 x 10 3 Kg/m3 

Then: - _ 18 (1.796 x 10~ 5 kg/m sec) = 2 57 x 1 

1.4 x 10 3 kg/m 3 (.3.x 10 6 m) 2 

For an estimate of the particle lag, u, where u = 
Ug-Up, consider the flow in the nozzle of the jet. As- 
suming the flow velocity increases from about 0 to 15 
m/s within .15 m then Ug may be approximated as a 
linear function along the nozzle, or. 

Ug = °=x where <= = 100 (5) 

Returning to (4) and making the approximations of 
a steady one dimensional flow with u << Ug then 
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ag = ug g + Ku (6) 

Substituting (5) into ( 6 ) 

9 . d u 

= z x = «x - 3 — + Ku 

ax — 

This can be integrated, and assuming no initial lag at 
x = 0 , then 

cc ^ c: 

u = x = Ug 

« t K « + K 


Thus at the exit plane of the nozzle, the lag can be 
estimated as: 


_ 1 x 10 2 Ug 

(1 x 10 2 ) + (2.57 x 10 6 ) 


4 x 10“ s Ug 


The particle lag is thus shown to depend on the 
velocity and viscosity of the flow, and the density and 
diameter of the particles. In this investigation, the 
lag is quite negligible. 

Next consider a simple turbulent flow whose 
correlation function is assumed to be 


R(t) 


u(t) ,u(t + t) 


u" 


exp (- 2 irf * t) 
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where f^ is a characteristic frequency of turbulence. 
The ratio of the rms particle velocity fluctuation/ 
0 .el, 2 ) 2 to the rms fluid velocity fluctuation (u 2 ) 2 is 


(up 2 ) 3 * _ 1 

(u 2 ) 55 1 .+ ft/fp 


with fp - — 
2 IT 


Recalling 

fp 


2.57 x 10 s 
2ir 


4 X 10 5 HZ 


Hence , 


the ratio 


(up 2 ) 8 ’ 

(u 2 ) 85 


is a function of the characteristic 


frequency of the turbulence but for large K, the 
quantity tends toward 1.00. 

Essentially, the response of the particles act 
similar to a low pass filter. The 3dB point of the rms 
response occurs at approximately 4 x 10 s Hz. 
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APPENDIX C 


Phase locked loop processor 

Description 

The basic function of the phase locked loop is to 
detect and track small differences in the frequency and 
phase between an input signal and a secondary reference 
signal. A simple diagram for the phase locked loop 
is shown in Figure 66(a) . The PLL is a negative feed- 
back loop consisting of a phase detector (PD) , a low 
pass filter (LPF) and a controlled oscillator whose 
frequency is either a linear function of a dc voltage 
or current. The PLL in this investigation used a current 
controlled oscillator (CCO) . 

As the loop tracks the incoming signal, the phase 
detector compares the frequency of the input to the out- 
put of the CCO. A change in the phase of the input is 
detected by the PD which produces an output proportional 
to the sum and difference in frequency and/or phase be- 
tween the input signal and the CCO output. The low pass 
filter then eliminates the sum component of the PD's 
output. An error signal is produced by the LPF whose 
magnitude is directly proportional to the phase difference 
determined at the phase detector. The error signal is 
then used to adjust the frequency and phase of the CCO in 
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Figure 66(b) Linear Model of Phase Locked-Loop 









such a direction as to establish phase synchronization 

between the two signals at the phase detector. 

Consider the following linear model for a PLL system, 

as shown in Figure 66(b) . 

$i and <f>'D represent phase angles associated with 

the input and CCO signals. F(S) represents current 

transfer function of the LPF and K, and K are conversior 

d o 

gains of the PD and CCO. Consider the input signal, 

Vi (t) and the CCO output, Vo (t) at the PD: 

Vi (t) = Vi sin Wit 
V Q ft) * V Q sin (w Q >t + ♦.) 

where Wi , w Q and <j> are the frequency and phase characte- 
ristics of interest. 

The PD ' s output current signal is: 

(sin w i t) (sin (w Q t + $)) 


V c (sin w t) (sin w t + 
cos (a-b) - h cos (a + b) 


i e (t) = K V.V 
e i x o 


where K is a constant, 
l 


Case I : 


w i = W o = w 


i 0 (t) = K V . 


Recalling sin a sin h ~ h 
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Then: 


K V. V 
1 x o 


(t) « 


{cos^ - cos (2 w t + $)) 


The low pass filter removes the high frequency part 
and so: 

K V. V 
2 i o 

i (t) = i d = -x cos $ 

6 4 


where K is another constant. 

Z 

If w^ and w q are perfectly synchronized to the free 
running frequency then i^ will equal zero which im- 
plies Q must be 90° . If w^ and w q are slightly out of 
phase then i^ is propo'rtional to the phase difference 
centered about the reference angle of 90°. If w^ now 
changes from w q then i^ will, adjust to a non-z'ero value 
in order to correct w Q so that frequency lock can be 

maintained with w. = w . § will be shifted some amount 

i o 

A 4> from 90° : 

Then i, becomes 
d 

i, - K ViVo cos { *o ± A<f>) - ± K V. V sin A a 

d ji T 2 1 ° 

In addition to being an error signal, i^ represents 

the demodulated output of an FM input applied as V_^ (t) . 

Case II: w. 4= w 

x T o 

This case corresponds to an unlocked conditon for 
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the loop. Hence the phase angle is meaningless and can 
be eliminated by appropriately choosing the time origin. 
The PD's output is then: 

V.V K = . 

i (t) = — ~ ' " 3 (cos (w, — • w ) t - cos ,w. + w ) t)- 
s £ x a x o ’ 

After filtering: 

i (t) = i. (t) = K V, V cos (w, - w )t 
e a 310 X o 

where K is a constant. 

3 

Notice that i^(t) represents the signal current 

which is supplied to the CCO. If the amptitude of Vj. 

is sufficiently large and w^ does not exceed the dynamic 

range of the CCO then lock-in can be established. That 

is, the CCO's output frequency will be shifted from 

the free running frequency, w^-' by some Aw. If w. is 

beyond the dynamic range, then w should be shifted 

o 

closer to the expected w^ frequency. 

Specifications 

The PLL used in this investigation had the cap- 
ability of determining the mean frequency (velocity) in 
the flow from 600kHz (3m/sec) to 4 MHz (20m/sec) . The 
useful frequency range over which the PLL can track the 
input signal (i.e. Lock Range) is ± 36% of the centered 
frequency. The frequency range over which the PLL can 
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acquire an input signal after lock has been lost (i.e. 
Capture Range) is set by the low pass filter’s characte- 
ristics and has a value of ± 28% of the mean or centered 
frequency . 

For the calibration curves of the PLL and the 
techniques used for the actual acquisition of information, 
please see Reference 53. 
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APPENDIX D 



Since it is extremely difficult to directly measure 
the velocity in a direction perpendicular to the direction 
of the mean flow, the following scheme was devised. 

Step 1.: Measure velocity at a positive (+) and 

negative (-) angle with respect to the 
velocity U at -a point in the flow field. 

For simplicity, = -0 « 45° 

Step 2:- For mean velocity: 


U+ = U cos0 + V sin 9 
1 1 



/2 

2 
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and 


Step 3 


then 


and 


V = (3+> ~ ( u ~) 

/2 

: For turbulent instensity: 

Xi+2 _ m >/ 2~ + v y ^~~ . 2 
'2 2 ’ 

U-2 - m _ v AAA) 1 

' 2 2 ' 

u+ 2 = {U2 + 2uv + V2) 

u- 2 = (u' 2 - 2uv + v 2 ) 

V 2 = (u+ 2 + u'- 2 ) — U 2 
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